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tion.  'Ihe  purpose  of  thv  present  research  is  to  isolate  the  effects  of  curvature 
(swirl)  on  the  turbulence  and  hence. the  transport  properties. 
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developing-  axial  and  decaying  tangential  velocity  fields  using  isothermal  air 
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ABSTRACT 


An  inport  ant  feature  of  many  flews  encountered  in  practice, 
(such  as  in  turbomachinery) ,  is  the  fact  that  the  streamlines  may  be 
curved,  thereby  introducing  a  pressure  gradient  in  the  direction 
perpendicular  to  the  main  flow  direction.  The  purpose  of  the  pre¬ 
sent  research  is  to  isolate  the  effects  of  curvature  (swirl)  on  the 
turbulence  and  hence  the  transport  properties. 

The  experimental  effort  is  concerned  with  mapping  out  the 
details  of  the  developing  axial  and  decaying  tangential  velocity 
fields  using  isothermal  air  as  the  working  fluid  in  an  annulus  with 
a  single  diameter  ratio  (d^/dg  =  0.4)  and  at  a  single  bulk  Reynolds 
number  (Re  =  1%/y  =  130,000.)  In  this  report,  interest  is  centered 
an  a  critical  discussion  of  the  data  reduction  techniques  for  ob¬ 
taining  the  radial  variations  of  the  axial  and  tangential  momentum 
diffusivities. 
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EXPERIMENTAL  'TURBULENT  VISCOSITIES  FOR. 

SWIRLING  FLOW  IN- A  SrATICNAKT  AMNULUS. 

by  C.  J,  Scott  and  D.  P„  Rask 
•PTOODUCTION 

-  An  important  feature  of  many  swirling  flows  encomtered  in  practice  (such 
as  in  trabomachiiiery)  is  the  existence  of  strongly  curved  streamlines  that  intro¬ 
duce  a  pressure  variation  in  a  direction  perpendicular  to  the  main  flow.  An 
engineer’s  principal  hope  of  tn&srstaiiding  such  in  obviously  complex  flow  system 
is  through  a  synthesis  of  some  simpler  modes  into  analytically  tractable  forms 
of  else  through  laboratory  simulations.  The  former  approach  has  only  been  fruitful 
in  laminar  swirl  situations  where  the  transporc  properties  are  known.  The  avail¬ 
ability  of  .digital  computers  and  refined  numerical  schemes  have  only  recently 
resulted  in  a  clearer  understanding  of  some  of  the  constant-property,  laminar, 
rotating  flow  systems.  On  the  other  hand,  rotating  turbulent  flows  are  diffi- 
■  cult  to  treat  because  of  the  lack  of  accurate  turbulent  transport  coefficients. 

For  the  relatively  simple  case  ,  of  a  turbulent  core  vortex,  the  values  of  turbu¬ 
lent  viscosity  suggested  in  the  literature  vary  ten-fold. 

The  experimental  effort  reported  here  is  concerned  with  mapping  out  the 
details  of  the  developing  axial  and  decaying  tangential  velocity  fields  using 
isothermal  air  as  the  working  fluid  in  an  annulus  with  a  "Single  diameter  ratio 

V''  — 

(dj/dQ  =  0.4)  and  at  a  single  bulk  Reynolds  number  (Re  =  Ut^/v  =  130,000).  In 
this  report,  interest  is  centered  on  a  critical  discussion  of  the  data  reduction 
techniques  for  obtaining  the  radial  variations  of  the  axial  and  tangential  momen¬ 
tum  diffusivities.  Thejannular  geometry  was  selected  because  the  inner  and  outer 
walls  produce  different  effects  on  the  swirl  turbulence. 

Consider  a  high,  axial  Reynolds-number,  flow  through  an  annulus.  If  both  walls 
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are  stationary,  then  if  a  wheel-like  rotation  (solid-body)  is  initially  imparted 
to  the  flow,  near  the  inner  wall  the  positive  radial  gradient  of  angular  momentum 
is  known  to  have  a.  stabilizing  effect.  Near  the  fixed  outer  boundary,  a  destabili 
ing  effect  is  encountered.  The  net  effect  of  the  centripetal -acceleration  field 
is  to  exert  a  stabilizing  influence  on  the  flow  which  depends  on  the  shape  of 
the  angular -momentum  distribution.  Local  stability  variations  modify  the  turbu¬ 
lent  structure  mainly  through  the  production  of  turbulent  energy  such  that  the  ;tur 
bulent  transport  coefficients  are  inhibited.  Near  the  outer  wall  -  the  unstable 
case  -  the  net  effect  of  the  centripetal  acceleration  field  is  to  increase  the  tur 
bulent  production  and  promote  turbulent  transport. 

The  characteristics  of  swirl  flows  have  been  investigated  during  the  last 
three  decades  and  are  still  the  subject  of  extensive  research.  Much  of  the  work 
was  simnarized  in  a  symposium  [1]*  on  "Concentrated  Vortex  Motion  in  Fluids,"  or¬ 
ganized  by  the  International  Union  for  Theoretical  and  Applied  Mechanics  and  held 
at  the  University  of  Michigan,  Ann  Arbor,  Michigan,  in  July,  1964  as  summarized 
in  J.  Fluid  Mech. ,  Vol.  21,  No.  1965.  Later  in  1966,  a  book,  Vol.  7  of 
"Progress  in  Aeronautical  Sciences,"  [2]  reports  in  more  detail  matters  discussed 
at  that  particular  synposium  and  presents  revised  versions  of  summary  reports  pre¬ 
sented  there. 

A  few  practical  applications  of  swirling  flows  are  listed  below: 

1)  Heat  transfer  per  unit  area  as  well  as  ratio  of  heat  transfer  to 
friction  loss  is  increased  if  rotation  is  superimposed  on  axial  flow 

4 

through  tubSs. 

2)  In  pipe  flow,  twisted  tapes  have  been  inserted  to  augment  heat  transfer. 
Applications  here  include  throat  regions  of  regeneratirely-cooled  rocket 
nozzles  (where  local  hot  spots  must  be  avoided)  or  in  high  performance 


•*  Numbers  in  brackets  designate  References  at  end  of  report. 
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•  heat  exchangers. 

3)  A  common  instance  .of  flow  between  concentric  cylinders  occurs  in  rotating 
machinery*  where  the  heat  transfer  characteristics  of  the  air  gap- are 
frequently  the  least  understood-  of  the.  many  heat  flow  paths  in  electric 

--  rotating  equipment. 

4)  Ip.  gaseous  nuclear  reactor  Tocket.  motors,  swirling  jtotion  achieves  a 
longer  hold-tp>  time  -of  the  fissionable  ;gas.. 

5)  In  .plasma :  gene ratois ,  a  swirling  motion  is  superimposed  on.  tit  axial 
gas  flow  through  the  nozzle  in  order  to  move  the  arc  attachment  point 
on  the  anode  circumferentially  and,  thereby,  avoid  overheating  the  anode 
material. 

6)  In  advanced  conbustion  chanbers,  swirling  motion  of  the  gases  improves 
flow  stability  and  combustion  efficiency  while  reducing  the  heat  loss 
through  the  walls. 

7)  With  increased  levels  of  tenperature,  cooling  must  be  applied  to 
rotating  and  stationary  conponents  of  gas  turbine  disks  and  shrouds. 

From  a  heat  transfer  standpoint,  this  poses  the  question  of  predicting 
the  heat  transfer  in  passages  enclosed  by  surfaces  which  are  partially 
rotating. 

8)  Swirling  flows  are  used  to  separate  particles  such  as  dust  separators 
in  the  air  intakes  oi  helicopters  hovering  over  dusty  fields. 

9)  Ranque-Hilsch  tube  cooling  devices  are  now  in  production.*  based  on  the 
energy  separation  of  a  gas  stream  into  a  portion  with  higher,  and 
another  portion  with  lower,  tenperatures  than  the  inlet  tenperature. 

10)  In  certain  flow  regimes,  the  effect  of  swirl  is  to  improve  diffuser 
performance;  while  in  others,  the  performance  is  reduced. 

11)  In  a  typical  submarine  condenser,  the  turbogenerator  exhaust  inlet  is 
located  on  top  and*  at  one  end  of  the  cond^tfer  but  is  offset  from  the 
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synmetfical  plane.  The  vaporenteringthe  condenser- through  the  inlet 
is-inalcgous-  -fovi  ’jit'  StaSaipg<.t4ngCTtially  antd-*an--®aulus.  Once 
within'  the  annulus  ,  the  vapor  toads  to  raove  spirally  toward  the  opr 
posite  end  of  the  cylinder.  -  ' 

12)  Che  mght  add  to  this  list  atmospheric  and'  geophysical  vortical- 

phenomena ’such  as 'dust  devils,  waterspouts,  tornadoes,  whirlwinds,  etc. 

Swirling,  flows  Juive  been  classified  by  King,  RotMus,  and  Kermpde  [3]  accord¬ 
ing,  to  the  boimfery  condition;  in  the  following  fashion: 

1)  Unconfined  swirling  flows. where  the  wall  effects  are  negligible; 

Exanples  are  swirling  free  jets,  rotary  gas  burners,  trailing  vortices 
from  aircraft  wings,  and  swirling  compressible  nozzle  flows. 

2)  Small  L/D^  swirling  flows  in  short,  large-diameter,  charters  where  side¬ 
wall  effects  strongly  interact  with  the  swirl  to  produce  significant 
secondary  flows.  An  example  here  is  core -flow  containment  in  gaseous 
nuclear  rockets. 

3)  Large  L/D^  swirling  flows  in  tubes  where  circumferential  wall  effects 
interact  strongly  with  the  swirl,  flow. 

Revolving  fluid  flows,  in- which  the  tangential  velocity  distribution  w(r) 
is  described  by  wr  =  constant,  are  normally  classified  as  vortex  flows.  The  revolv¬ 
ing  fluid  is  caused  to  move  inward  toward  the  axis,  of  revolution.  If  the  fluid 
friction  is  low,,  or  if  the  radial  flux  of  angular  momentum  is  high,  the  fluid  exer. 
cutes  a  vortex  motion.  These  flows  are  normally  in  contact  with  one  or  more  solid 
boundaries.  Flows  may  be  rotating, over  stationary  walls  or  inside  confining 
stationary  walls.  The  boundaries  may  be  rotating  and  inpart  swirl  to  the  main 
body  of  fluid  by  turbulent  diffusion.  In  the  case  of  a  grounded  storm  funnel, 
both  the  fluid'and  the  surroundings  are  rotating. 

For  swirling  internal  flows  in  long  stationary  ducts,  the  axial  flow  (u) 
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assumes  a  fully  developed' velocity  distribution  at  large  distances  downstream- and 
the  swirl  (w)  decays  asymptotically  to  zero.  If  wail  rotation  is  present *  a  fully 
developed  swirl  profile  is  produced.  If  the  swirl  ratio  is  sufficiently  large*, 
stagnant  regions  arid  regions  of  reversed  flow  may  exist  in  the  entrance  section. 
This  phenomenon,  known  as  vortex  breakdown ,  niay  occur  in  swirl  diffusers  and 
influences  their  performance.  It  helps  to  stabilize  flames  arid  alters  the  flow 
frrai  rotating  jets.  It  Has  been  proposed  to  contain  a  volume  of  fluid  within  a 
body  of  different  fluid  with  relatively  little  mixing. 

Bossel  [4]  investigated  laminar  vortex  breakdown  flows.  One  special  case 
involves  swirling  flow,  in -a  cylindrical'  stream  tube.  These  solutions  are  char¬ 
acterized  by  a  maximum  swirl  angle  of  &  ~  62,5°.  At  larger  swirl  angles,  a 
stagnation  point  and  a  reversed  flow  exist  along  the  axis  of  the'  tube.  The  pres- 
erit  experimental  studies  are  confined  to  the  subcritical  case  where  vortex  break¬ 
down  does  not  occur  because  of  the  difficulty  in  deducing  distributions  of  turbu- 
lent  transport  properties  with  reversed  local  axial  flows. 

Another  interesting  feature  of  swirling  flows  is  related  to  their  hydro¬ 
dynamic  stability.  Since  particles  in  revolving  motion  tend  to  conserve  their 
angular  momentum,  the  stability  criterion  first  proposed  by  Rayleigh  [5]  suggests 
that  flows  with  positive  radial  angular  momentum  gradients  are  stable  while  flows 
with  negative  gradients  are  unstable.  Rayleigh  describes  the  equilibrium  behavior 
of  a  fluid  particle  near  the  outer  wall  in  a  curved  channel.  Upon  being  displaced 
to  a  larger  radius,  the  fluid  possesses  a  larger  angular  momentum  than  its 
neighbors.  The  centrifugal  force  on  the  displaced  particles  (-pw  /r)  will  be 

'  ■  H  ■'  1  *  1  1  '  « 

*  The  term  swirl  rate  (ratio)  used  here  is  subjective  in  nature,  (hie  could  say 
a  swirl  flow. is  of  a  high  rate  if  in  the  inlet  region  -(u^  and 

are  the  maximum  axial  and  tangential  velocities).  The  largest  swirl  rate;  reported 
here  has  a  ratio  of  1.6  at  the  first  data  station  (x/r^  =  1.7). 


greater  than  the  centripetal  pressure  gradient  existing  at  the  new  location  and 
the: pa.  tele  will  tend  to  be  displaced. to  a  still  larger  radius.,  :Ihere£6re,  a 
radial'  displacement  is  unstable  near  the- outer  wall.  Near  the  inner  wall,  the 
converse  is  tru$.  _  Rp  tting  flows  with  solid  body  rotation  (w*  =  .cr)  are-  very 

't  .  \  '  • 

stable,,  "free"  vortex  flows  (wr  =  c)  are  neutrally  stable,  and  the  boundary  layer 

v .  .  i 

region- inside  the  channel  wall  is  extremely  unstable.' 

In  1932  and  1936,  Taylor  [6,  7;;.-81  investigated  fully- developed  rotating 

turbulent  flow  between  concentric  rotating  cylinders.  In  Taylor's  rotating  cylin- 

!  '  1 

der,  the  Taylor  vortices  appeared  when  the  outer  cylinder  was  stationary  and  the 
inner  cylinder  was  rotating.,  In  that  case,  the  angular  momentum  distribution  is 
unstable.  In  Gortler’s  experiments  [9]  with  flews  over  concave  and  convex  walls, 
the  3-D  disturbances  appeared  only  for  flow  over  concave  walls.  \  ■ 

In  1935,,  Wattendorf  [10]  attenpted  to  isolate  the  effect  of  streamline 
curvature  on  turbulent  flow  by  using  a  curved  channel  of  constant  curvature  and 
cross  section.  He  found  that  for  fully-developed,  curved,  internal  flows,  the 
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turbulent  viscosity  data  were  less  than  that  for  straight  flow  near  the  inner 

’  1  \  .  \ 

(convex)  wall  and  greater  near  the  outer  wall  fin  accordance  with  Rayleigh's 
criterion  of  stability.  Viscosities  near  the  outer  wall -were  up  to  four  times 
corresponding  inner  wall  values.  Eskinazi  and  Yeh  [11],  and  later  Margolis  and 
Lumley.  [12] ,  demonstrated  that  the  -turbulent  structure  is  dependent  on  the  shape 

of  the  angular  momentum  profile.  Their  results  with  curved  channel  flows  led  \ 

*  5 

to  the  conclusions  that  on  the  convex  wall  the  ^gradient  of  angular  momentum  is 

\  ! 

positive  and  turbulence  is  suppressed;  whereas  on  the  concave  wall,  the  gradient  - 

of  angular  momentum  is  negative  and  turbulence  is  promoted.  In  1968,  Bradshaw  [13] 

concluded  that  the  effect  of  streamline  curvature  on  the  turbulence  is  ^probably 
’  ,  .  1  1 
more  important  than  the  effect  of  streamline  curvature  on  the  mean  equations.  To 

conclude,  the  chain. of  reasoning  is  that  streamline  curvature  influences  the  flow 
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^stability,  thereby  modifying  the- local  turbulence  properties.  Swirling  flow 
in  an  annulus  is  adopted  here  as  the  experimental  geometry  because  it  provides 

variable  streamline  curvature  and  because  annular  flows  apjlear  to  develop  and 

—  *  _  .  *. 

dec ay  more  rapidly  than. pipe  flows. 

Four  pertinent  survey  articles  on  swirling  flow  are  by  Westley  [14], 

Rrsith  '[15]  >.  L&vin  and  Fejer  [16],  and  Gairbill  and  Bundy  [17],  In  1937,  G.  J. 
Ranque  [18]-  applied;  for  a  vortex  tube  patent.  .In  1944,  Hilsch  developed  a  vortex 
cooling  de\rice  wHdi  freqiKntly  .bears  his  name.  In  1946-48, '(Milton  [19],  and 

'  ..  .  .  •  j 

later, -Reed;  [20], tpointed.out  the!  important  phenomena  occurring  in  vortex  tubes. 

”  ,*  '  '  '  \  - 

In  ,1948,  Kasssier  and' Ifeo^mschild  [21];  explained  the  energy  separatipn  assuming 

that'  the^pptential  vortex  (wr  »  const)  initially  formed  at  the  nozzle  inlet  is 
transformed  into  solid  body  rotation  (w/r  =  const)  as  the  flow  proceeds  axially 
toward  the  outlet,  .A'.sirjilax  analysis  was  given  by  Deissler  and  Perlmntter  [22,23] 
who  proposed-- that  the  axial  flow  Was  necessary-  in  the  energy  separation  process 
because  most  or  the  fenperature  separation  . occurs  where  tire  bulk  flew  changes  from 
radial-  to  \axial.  \ 

I 

.  Jn  papers  by  Dsissler  and  Fcrlmutter  [22]  and-'  Donalds.on  and  Sullivan  [24], 
the  flow  is  assumed  to  be  twd-diiipnsianal  or  have  ah  axial  velocity  proportional 
to  the  axial,  coordinate.  lift  -11113  -case,  Vthe\shaps  .of  the  tangential  velocity 
distribution  is.  controlled -by.  the  radial.  Reynolds  ntsiber.,  -Re  »  and  is,  independent 
of  the  fahgentiai  velocity  magnitude,  For  very  lew  Re  ,  the  tangential  velocity 


is  found  to  hew  =  r  (solid-body  rotation))  As  .Rer  increases,  the  velocity 
profile  approaches  w  c  r- 1  (free  vortex  flow)  except  nfear  the  axis  where  solid 
body  rotation  always  occurs;  Kreifh  and  Margoiis  [25],  Snithberg  and  Landis  [26], 
and  Thorsen  and  Landis  [27]  studied  the  effects  of  turbulent  swirling  motion  cn 
heat  transfer  from  tubes.  Krassner  and  Knoerschild  [21],  Hartnett  and  Eckert  [28], 

I 

Deissler  and  Perlmqtter  [23],  Ragsdale  [29],  Keyes  [30],  Kerrebrock  and  Meghreblian 


i 


8, 


[31],  Reynolds  [32] ,,  Kendall  [33],  and  SibiiCcin  [34],  all  studied  forms  of  the 
Rahque.-Hilsch  tube. 

Collatz  and-Gortler  [35]  and  Talbot  [36]  analyzed  swirling  flow  in  pipes 
•by  linearizing  the  tangential  momentum  conservation  equation  and  employing  a 
constant  turbulent  viscosity.  The  simplified  equation  is  often  referred  to  the 
"swirl"  equation.  Kreith  and  Sonju  [37] ,  Rochino  and  Lavan  [38] ,  and  Scott  [39] 
all  extended  the  method  by  varying  initial  conditions,  using  variable  properties 
and  annular  geometries.  Muslof  [40],  measured  the  decay  of  turbulent  swirl  in  a 
stationary  pipe,  Persen  [41]  studied  swirling  inlet  flow  in  a  tube  using  a 
boundary  layer  approach.  Rask  and  Scott  [42]  studied  the  decay  of  turbulent  swirl 
in  an  annular  duct.  Boemer  [43]  employed  the  time-averaged  Navier  Stokes  equa¬ 
tions  for  turbulent  flow  for  the  annular  geometry.  In  Boemer's  work,  the  effects 
of  turbulent  transport  are  incorporated  through  the  use  of  two  apparent  turbulent 
viscosities  which  are  computed  separately  for  use  in  the  axial  and  tangential 
momentum  equations.  Best  results,  as  measured  by  agreement  with  the  data  of 
Rask  and  Scott  [42] ,  were  obtained  with  a  tangential  viscosity  model  that  links 
the  viscosity  to  tangential  velocity  via  the  average  flow  angle,  the  wall  shear 
velocity,  and  the  distance  cf  the  zero  shear  location  to  the  wall. 

From. a  review  of  the  literature,  it  is  evident  that  only  a,  limited  number  of 
previous  investigations,  either  experimental  or  analytical,  have  been  devoted 
to  vortex  flows  in  concentric  annuli.  A  fair  number  of  investigators  have  studied 
the  flow  between  concentric  rotating  cylinders  in  the  .absence  of  an, axial  flow, 
but  this  is  relatively  far  removed  from- the  present  investigation. 

Palanek  [44]  investigated  swirling  flow  in  au  annular  duct  of  diameter  ratio 
d^/dQ  =  0.307.  Palenek's  study  was  initiated  as  an  extension  of  an  earlier  work 
by  Taicott  [45]  on  essentially  the  same  apparatus.  Talcott  investigated  heat 
transfer  characteristics  in  straight  flow  and  also  briefly  in  swirling  flow.  The 


swi$ 1  used  in-Palenk’s  study  were  so  extreme  that  reversed  flow  is  observed 

iii.eadi;4j3£e.'.  _ 

As.t^ly  cicely. related  to  tiurt  reported  here  was * done  by  Yeh  [46] ;  The 
uaajor  differences ;^x©  that. his  test  sectiqh-was  of  larger  hydraulic  .diameter 
and  diameter  -ratio,;  arid 'was  relatively  shorter  "than  the  present,  test,  section  In 
terms  of  hydraulic  diameters.  In  an  annulus,  the  hydraulic  diameter  is  given  by 
=  dQ  -  d:..  The  pertinent  dimensions  of  Yeh’s  apparatus  were  -dj/d  =  0,5  and 
=  5  inches.  Yeh  reported  lbcal  skin  friction  coefficients  at  the  inner  and 
outer  walls  but  the  investigation  is  somewhat  inconplete  in  that  only  straight 
flow- and  one  swirl  rate  were  investigated.  Acharyn  and  associates  [47]  also 
analyzed  the  data  that  was  obtained  by  Yeh.  Hie  shear  stress  results  are  in  poor 
agreement  with  those  reported  by  'Yeh. 


Turbulent  Transport  Properties 

For  laminar  flow  along  a  curved  path,  x£  =  p(3u/3r  -  u/r).  It  is  a  custo¬ 
mary  but  questionable  practice  to  use  an  analogous  expression  for  the  turbulent 
shear  stress,  i.e.,  t^.  =  pet(3u/3r  -  u/r),  where  et  represents  ah  apparent  tur¬ 
bulent  coefficient  of  kinematic  viscosity  or  the  "exchange  coefficient"  of  the 
mixing  process.  Using  this  approach  xt  «  0  when  (3u/3r  -  u/r)  =  0.  Kline  [13]' 
has  pointed-  out  that  if  an  experimental  shear  stress  is  divided  by  an  experi¬ 
mentally  determined  gradient  of  velocity,  the  ratio  can  go  to  plus  or  minus 

infinity.  This  leads  to  difficulties  in  plotting  experimentally  determined 

\  * 

transport  coefficients.  For  rectilinear  flows,  Prandtlts  mixing  length  and 
Taylor's  vorticity  length  hypotheses  give  x^  =  i.2|3u/3y|3u/3y.  Again  the  shear- 
stress  depends  on  the  velocity  gradient,  it2  is  the  square  of  the ^mixing  length 
and  a  function  of  the  space  position.  Von  Kantian  tried  to  relate  turbulent  trans¬ 
port  properties  to  local  conditions.  His  similarity  of  turbulent  motion  concept 
yields  xt  =  const  (3u/3y) V(32u/3y2)2  for  rectilinear  flows.  This  is  also  a 


jpradientmodel.  * 

Pranitl  [48]  ^plied  Rayleigh ' s  instability  reasoning  to  developed  curved 
turbulent  flows-,  assuming  the  local-  angular  momentum  (ur)  constant  during  dis¬ 
placement  of  fluid ^  particles  -perpendicular  to  mean  streamlines  and -leads  to 
;|rurj'?2.9(ru)/3r 

xr  »  .  c£v>f  -  £  GO 

viiere  (3»V3r  +  u/rj  represents  the  vorticity  at  thepqint,  v  the  normal  velocity 
and' £  the  mixing  length.  This  formula  is  identical  to  Prandtl’s  formula  for  shear 
stress  in  parallel  motion-except  the  velocity  slope  is  replaced  . by  (9u/3r  +  u/r)  . 
Equation  (1)  leads  to- a  vanishing  turbulent  shear  stress  for  irrotational  mean 
flows --a  behavior  contrary  to  the  experimental  observations  of  Taylor,  Wattendorf , 
and  Rask  and.  Scott. 

Kinney  [49],  in  dealing  with  plane  rotating  turbulent  flows,  in  the  absence 
of  ah  axial  component,  applied  an  extension  of  von  Karman’s  similarity  hypothesis. 

A  family  of  similar  velocity  profiles  was  generated  which  includes  both  the  rota¬ 
tional  solid  body  case  and  the  irrotational  free  vortex  case.  A  universal  constant 
appearing  in  the  equations  was  evaluated  using  G.  I.  Taylor's  rotating  cylinder 
velocity  and  wall  shear  stress  measurements.  Kihney  relates  the  local  angular 
velocity  of  the  fluctuating  flow-  to  the  local  angular  velocity  of  the  mean  flow 
with  the  relation  |w'7rj  *  l  3(w/r)/dr,  and  concludes  that  in  plane-curved  tur¬ 
bulent  flows  which  possess  universal  similarity,  the  angular  velocity  is  a 
transferable  "quantity,  just  as  linear  momentum  is  a  transferable  quantity  in 
rectilinear  flows. 

Rochino  and  Lavan  [38]  added  the  axial  dimension  to  the  problem  in  investi¬ 
gating  turbulent  swirling  flow  in  stationary  ducts  by  using  both  Taylor's  modified 
vorticity  theory  and  von  Kamrn's  similarity  hypothesis  [SO].  Three  possible 


It  is. not  yec  established  which  relation  is  Valid:  The  nqiinal  procedure  seems  to 
adopt  an  analogy  to  the  laminar  equation  vr’w1  =  pet(3w/3r  -  w/r)  where  theeddy 
di£rusivity  for  momentun  is  related  to  the  mixing  length  £.  by  =  const  zz  ’ 

(3w/3r  -  w/r).  The  radial  dependence  of  the  mixing  length  is  then  the  focal  point 
of  the. controversy,  hi  the -core  flow  away  ;from  any  wall  influence,  z  =  Jer, 
i.e.,.  the  mixing  length  increases  at  a  rate  proportional  to'  the  distance  from  the 
axis  of  -symmetry.  Near  the  walls,  the  mixing  length  increases  at  a  rate  propor- 
tional  to  the.  distance  from  the  wall.- according  to  Prandtl.  Since  the  annular 
geometry  has  two  walls  with  opposite  wall  behaviors,  it  provides  a  versatile  test 
geometry  to  study  the  anomalous  shear  stress  assumptions. 

Both  swirl  and  ah  annular  geometry  generate  asymmetric  mean  axial  velocity 
distributions  such  that  the  locus  of  zero  axial  shear  and  the  point  of  zero 
velocity  gradient  do  not  coincide.  The  computed  values  are  negative  in  a  small 
range  of  radial  distances,  Both  Wattendorf  [10]  and  Eskinazi  and  Yeh  [11]  showed 
that  in  a  fully-developed  curved  channel  having  an  unsynwetrical  profile,,  the 
location  of  the  zero  turbulent  stress  (-pu’v1)  did  not  occur  at  the  .point  where 
the  velocity  derivative  (3u/3r)  was  zero  nor  where  the  viscous  shear  was  zero 
[u(-3iV3r  -  u/r)  =  0]. 

The  difficulties  of  having  opposite  signs  for  the  shear  stress  and  the 
velocity  gradient  are  discussed  at  great  length  in  [13].  Alternate  proposals 
are  given.  Of  the  three  shear  expressions :  t  =  pei  3u/3y,  r  =  pe2/lfr7  3u/3y, 
x  =  pe3q2,  (q2  =  uTuT) ;  the  third  model  proposed  by  Townsend,  Bradshaw,  and 
McDonald  [13]  is  the  only  non-gradient  model.  This  model  requires  that  if  a  tur¬ 
bulent  kinetic  energy  exists,  a  shear  stress  exists.  Therefore,  the  production 


of  turbulent  kinetic  energy:  aliases,  increased  importance.  In  cylindrical 
x»  r»  0»  u>  v,  w,  coordinates,  the  production  of  total  turbulent  energy  is  [50] 
written  for  fully-develcped>.  isotropic,  swirling  flow  as 

.ot-|h]  (3) 

Ihe  production  is  positi\  o\er  the  entire  duct.  f  w  =  0,  or  with  solid 
body  rotation  w  =  cr,  the  normal.  xial  flow  production  results.  For  a  free  vortex 
w  “  c/f*  ?he  terra  3  (w/r)/3r  is  nev  iive  and  the  possibility  of  subtraction,  occurs 
depending  on  the  correlation  v'w'.  therefore,  swirling  flows  also  have  the  dif¬ 
ficulties  of  negative  "apparent"  tangential  turbulent  viscosities  defined^as 

_  v*w!  , 

6t  8W  _  W  W 

8r '  ”  ,r  . 

It  is  informative  to  examine  the  equations  used  to  deduce  the  turbulent  shear- 
stress  and  eddy  diffusivities  when  developing  profiles  occur. 

The  Governing  Equations 

Consider  a  cylindrical  coordinate  system  with  x,  r ,  and  e-  as  the  axial, 
radial,  and  azimuthal  coordinates.  Let:u,  v,  and  w  be  the  time  mean  velocities 
in  the  x,  r,  and  8  directions  and  u' ,  v,:,  and  w*  be  the  corresponding  turbulent 
fluctuation  velocities.  Writing  the  continuity  and  Navier-Stokes  equations  in 
this  notation,  taking  mean,  values  with  respect  to  time,  requiring  the  mean  motion 
to  be  steady,  the  mean  density  to  be  constant,  density  fluctuations  to  be  negli¬ 
gible,  and  symmetry  with  respect  to  e,  results  in  the  following  governing 
equations  [52]: 

Cohservation^of ‘ mass 
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Conservaiicn- of- 


axial  momenta: 


Conservation  of  Tadial  roraentta 
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Conservation  of  angular  .momentum 


fv'wn 


vlut  1 


The  resulting  set  of  turbulent  Navier-Stokes  equations  is  non-linear  and; 
indeterminate  due  to  the  turbulent  shear  terms  at  the  extreme  right  in  equations 
6-8.  Since  there  are  more  dependent  variables  than  governing  equations,  no 
general  method  of  solution  exists.  We  shall  now  form  integrals  of  the  equations 
of  motion  in  the  following  manner. 

(a)  Equation  6  is  multiplied  by  r  and  integrated  across  the  flow 
from  r  =  r^  at  the  inner  wall  to  r  =  r.  The  result  is 


where 


9u  — , r 

**3r  "  PU*F 


=  p(v  +  e 

rx-'sr 


is  the  axial  shear  stress.  The  tern  y32u/3x2  was  omitted  in  equation  9  in 
accordance  with  an  order  of  magnitude  analysis  using  the  data  taken  in  this 
study  which  revealed  that  except  for  regions  close  to  the  walls  (closer  than 
the  closest  .experimental  point),  the  turbulent  viscous  terns  in  the  equations  of 
motion  are  on  the  order  of  SO  times  the  molecular  terms.  Coupling  this  fact 


with  the  normal  assuqitian  that  in  the  sidewall  boundary  layers  (d/.dx  «  3/3r), 

V2  2  sVSr2  +  l/r-3/3r  applies  accurately  for  molecular  quantities.  Equation  9 
is  made- non-dimensional  by  dividing  by  1/2  pu2  r  (u  is  the  average  axial  velocity) 


IX 
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(9*) 


The  five  terms  on  the  right-hand  side  of  equation  (9')  represent 
.1)  The  axial  shear  on  the  inner  wall  of  the  annulus 

2)  Mean  velocity  changes 

3) .  Kinetic  energy  changes 

4)  Pressure  changes 

5)  Turbulence  intensity  changes 


Integrating  across  the  entire  annulus  yields  the  following  relevant  integral 
relationship 


r°v0 '  riX 


d 

3E 


(P  +  uz  +  u72)  rdr 


ri 


(10) 


where  the  left  side  corresponds  to  the  total  shear  force  per  unit  length  (t  <  0) 

jL 

and  the  right  side  represents  the  change  in  axial  flux  of  linear  momentum.  The 
turbulence  level  and  the  radial  pressure  variations  are  important  here  because 
of  the  tangential  velocity. 


15 


-  "  @0-  Equation  8  is  nultipliedbyr2  and  integrated  from  the  inner  wall 

outward  (agshunsglecting  pa^/Sx2) 
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Hie  four  terns  on  the  right-hand  side  of  equatian.  il  represent 

1)  The  tprque  per  wit  length  on  the  inner  wall  of  the  annulus 

2)  Change  in  axial  flux  of  angular  momentum  between  inner  wall  and 
radial  location  r 

3)  Effects  of  axial  acceleration  on  the  axial  flux  of  angular  momentum 

4)  Hie  turbulence  shear  stress  rQx  =  u'w' 

Integration  over  the  entire  annular  gap  yields 

fT° 

!  r , 
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(iw  +  u'w')r2dr 
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(12) 


Hie  left-hand  side  is  the  total  torque  per  wit  length  while  the  right- 
hand  side  is  the  change  in  the  axial  flux  of  angular  momentum.  Equation  11  is 
made  non-dimensional  by  division  by  1/2  pTr2 
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(11') 


yiPr2  3x 


uTwrr2dr 


Ihese  equations  do  not  assume  that  the  direction  of  ' the  shearstress  is  the 
saae  as  the  resultant  velocity  or  resultant  velocity-gradient. 

Ihe  radial  aonentun  equation  integrates  directlyto 
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(13) 


According  to  Hinze  [51]  ,  the  tern  aCu^/Sx  is  of  smaller  order  of 
magnitude  than  the  other  -turbulence  terms. 

Finally,,  the  integral  expression,  for  mass  conservation  results  from 
direct  integration  of  its  differential  fern  (equation  5). 


ZV 


d 

clx 


undr 


(14) 
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No  turbulence  measurements  were  made  during  the -present  investigation. 
Equations  9’  and  11'  were  used  (neglecting  the  fifth  terra  in  equation  9*  and 
the  fourth  term  of  equation  11')  to  obtain  Tadial  profile.?  of'  axiai;  and'  tangential 
shear.  Axial  wall  shear  values  were  obtained  using  a  Preston,  itbe  or  a  Ciauser 
chart.  Tangential  wall  shear  values  were  obtained  using  the  flow  angle  at  the 
wall,  i.e.,  xr0  =  tan  <j>. 
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Eddy  di^usivity  ncitels  have  proven  to  he  a  useful  tool— parti raiiariy>  in 
adapting  laminar  flow -  solution  techniques;  to  permit  .calculations  in  cqnplex 
tyrtnileiitr  flows.  Since  all  diffusivity  models  ultimately  depend  on  empirical 
methods;  their  valise  must  be  tested  by  coup arison. with  reliable  erqjeriifints. 

It  is  the;  intention' of  the  present  paper  to  describe  sidt^bie'  e^eriinents 
carried  out  at  the  Heat  Transfer  Laboratory  of  the  L^iiversity  of  Minnesota. 

The  present  eipeinmental  program  was  carried  out  in  an  open  circuit  wind 
tumel  designed  fir  swirl  flews.  A  detailed  description  of  the  facility  is  given 
byRask  and  Scott  in  [42].  A  schematic  drawing  is  shown  as  Fig.  1.  Pertinent 
photographs  of  the  experimental  equipment  are  presented  in  Figs.  2-5.  Initially 
room,  air  enters  the  inlet  radially  ’where  airfoil-shaped  guide  vanes  impart  a 
swirling  motion  to  the  flow  prior  to  entry  into  the  annular  test  section.  Upon 
emergence  from  the  test  section,  the  mean  velocity  head  is  recovered  in  a  conical 
diffuser.  The  swirl  is  then  removed  from  the  flew  in  a  cylindrical  duct  filled 
with  straightening  tubes.  The  flew  is  then  monitored  at  an  orifice,  controlled 
by  means  of  a  low  resistance  by-pass  and  finplly  flows  into  the  inlet  of  a  blower 
which  exhausts  to  the  atmosphere  outside  of  the  laboratory. 

The  inlet  section  receives  the  initial  radial  inflow,  imparts  a  swirling 
motion  to  it  and  finally  converts  the  flow  into  the  axial  direction.  The  inflow 
enters  between  two  converging,  saucer-shaped,  guide  plates,  each  48  inches  in  dia¬ 
meter  with  a  radius  of,  curvature  of  42  inches.  Twelve  airfoil-shaped  inlet  guide 
vanes  are  evenly  spaced  on  an  18-inch'  circimference.  The  guide  vanes  are  col¬ 
lectively  rotated,  thereby  producing  variable  swirl  velocities.  Tire  vane  angle 
setting,  corresponding  to  a  straight  axial  flow  is  designated  as  zero  degrees. 
Increasing  the  vane  angle  of  attack  increases  the  downwash  velocities  behind  a 
blade  and  increases  the  resulting  tangential  velocities.  The  airstream  acquires 
angular  momentum  in  passing  through  the  cascade  of  vanes  at  the  inlet  of  the 
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'  turmel>  -  If  no  ietardihg'  forces  are- present ,j  'the  moment  of  moirehtum.is'  constant 
-  .as  the  £tow  spirals  inward  arid  a- free  yb^ex  type*  of  flow-  is 'generated;  The 
actual  situation -is  complicated  by-the  drag^of  the  vanes  i  .the.  presence  of  wakes 
\  ^behind  the.  vanes  ,  the  sidewall  plate  drag,  arid  the- drag' along  the  iraier  wall 

I  .........  .  '  '  '  ..  '  \- 

of  theranhiilus  within  the  inlet.  .  ,  5 

As  the  flow<encounters  the -90-degree-bend^which.  converts  the  principal  notion 
into  ari  aiaal-  motion,  it  encounters  three,  ring-shaped-,  boundary  layer  trips 
(Figure  lc);.  'The  trips  serve  to  fix  the  point]  of' transition  from  laminar  ^to 
.  turbulent  flow  at  th^  entrance  -to  the  annulus,  arid  thereby  achieve-  a  steadier 
gross  flow.  Olson  and‘:Sparrow  fS2]-  and'Okiishi  and^Serovy  =[53]  found  that  in 
siiooth, rounded-entrance  annuli,  the.  loots  of  transition- is  asymmetric  and  unsteady. 

.  '■  .  \  .  .  ;  "  t  .  ■  : 

!fhe  irregular;axial  drift  of  the  transition  point  causes i an  unsteadiness  of  .the 
flow  which  was  easily  measurable  when  the  trips  are  relieved.  >. 

The  test  section- is  a.conceiitric,  -constant-area,  annular  duct*  ihe-outer 
tube  is  of  clear  plexiglass  with-. a  5. 0-inch  I.D.  and-  1/4-ihch  w^ll  thickness. 


A  stainless  steel  irlner  tube,  2. 00-inch  O.D.  and  1/16-inch  wall  thickness  is 

used.  The  inner  tube  can  be  translated  axially  or  removed  completely.  The 

\ 

overall  length  of  the  test  section  is  118' indies  of  39.3  hydraulic  diameters, 
l  ' 

The  inner  stainless  steel  tube, is  instrumented  along  a  single  ray  with  29 
static  pressure  taps.  Variable  flow  spiraling  prevents  repeated  roughness 
effects.  Six  probe  wells,.  Figure  lb, rare  located  in  the  inner  tube  in  order  to 
accept  the  cylindrical  pressure  survey  probe.  The  probe  wells  are  installed  90 

i  ’ 

degrees  from  the  static  .pressure  taps  and  at  longitudinal  locations  midway  between 
two  pressure  taps.  ^ 

The  conical  diffuser  commences  with  the  S-'inch  diameter  plexiglass  tube  and 
increases  to  'an  18.4-inch  circular  cross  section.  The  walls  form  a  total 
included  angle  of  7  degrees.'  Straightening,  tubes  made)  from  mailing  tube'i-inch 
in  diameter  and-18  inches  long  fill  the  diffuser  cross  section  upstream  of  the 


I 


1 


.22 


r*;. 
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meteri^oiific«.  Their  primary  . function -is  -tp- remoy^  the  swirl  si^.  .that  the 
'tpt^/duDmlr.-fLaw-.'cini-'be  -accurately  measured.  The  orifice, -was  calibrated  by' 
-numerically integrating- the.  axial  velocity; profiles.  The  -calibration  was  found 
to  , be  independent,  of  the>  magnitude  of  the  s^flr-demonstfating  the  effectiveness 
of  the'  diffuser-straightening-tiibe-bank  combination.  The.  mass  flow  rate  is 
controlled  by  :means  of  a  low- resistance  by-pass .  A  section  of  the  tunnel,  loca- 
ted:be.tween  the  orifice  aid- the  blower,  has  bleed  holes  in  it  idiith-can.be 
exposed  by-  means  of  a  rotatable  sleeve.  By,  .opening  the  -bleed  holes,,,  room  air 
caii.be  induced  into  the  system  downstream- of  the  orifice!- -thereby  reducing  the 
rate;  of  mass  (flow  through  the  test  section. 

ihe  blower  used  to  draw  air  through  the  tunnel  is  powered  hy,  a  3  h.p. 
motor  equipped- with,  an  adjustable  v-belt  drive.  The  blower  is  isolatedifrpm 
its  angle-iron. mounting  frame  by  four  .coil  springs  and  isolated  from  the, 
annular  channel  by  means  of  flexible  canvas  sections. 

The  cylindrical  pressure- probe  and  the  hot  film  probe  used  to  sense  the 
local  velocities  are  each  held  in  identical  brass  blocks  which .caq.be,  .displaced 
vertically  along  a  dovetail  track  by  a  micrometer  which  is,,  in  turn,  secured  to 
the  base  section  cx  .a  surveyor's  theodolite.  Either  probe- can  be  rotated' a  full 
360  degrees  and  the  rotation  angle  determined  to  within  five,  minutes.  The 
micrometer  has  a  2-inch  travel  and  is  graduated  in  O.OO.l-inch  increments. 
Comparisons  between  Pitot  cylinder  probe  and  hot  film  probe  performance  are  made 
in -Appendix  A.  , 


Measurement  of  the  Mean  Values  of  the  Local  Static  arid  Total  Pressures .  The 
radial  velocity  distribution  v  may  be  calculated  from  the  continuity  equation  5. 
This  relation,  written  for  rotationally  symmetric  flow  with  steady  mean  values, 
holds -without  deviating  more  than  a  few  percent  from  the  complete  turbulent 
equation  of  continuity;  Integrating  equation  5  yields 


rv 


■r 

|-(ur)dr 


(IS) 


The  radial  velocity  is  identically  zero  at  both  the  inner  and  the  outer  walls. 

Only  near  the  entrance  are  axial  variations  of  .the  axial  velocity  of  signifi¬ 
cance.  The  average  axial  velocity  remains,  constant.  Therefore,  the  radial  com¬ 
ponent  of  the  velocity  is  everywhere  small  and  the- direction  of  the  flow  can  be 
described  by  a  single  aiigle  $  between  the  total  velocity  and  the  axial  direction. 

In  surveying  the  pressure  fields  in  internal  flows  with  stperinposed  secondary 
flows,  most  investigators  have  observed  that  the  flow  is  sensitive  to  distur¬ 
bances  resulting  from  the  presence  of  a  conventional  Pitot  probe. +  These  dis¬ 
turbances  are  minimized  when  a  Pitot  cylinder  is  used  which  is  oriented  along  a 
radius  and  is  of  sufficient  length  to  span  the  cross  section.  Ihe  principal 
design  and  construction  problem  of  this  type  of  probe  is  in  establishing  a  known 
angle  between  the  static. holes.  The  present  Pitot  cylinder  (Fig.  6)  was  made 
from  0.083-inch  O.D.,  0.063-inch  I.D.  hypodermic  tubing.  In  order  to  use  the 
probe  as  a  flow  direction  sensor,  two  pressure  taps,  0.020-inch  in  diameter,  were 
located  90°  apart  on  the  circumference.  By  locating  the  pressure  taps  45°  on 
either  side  of  the  stagnation  point  (Fig.  7) ,  maximum  sensitivity  to  the  flow  di¬ 
rection  was  obtained.  Thom  [54]  demonstrated  that  the  wall  pressure  sensed  in  a 
hole  drilled  in  a  cylinder  is  not  the  pressure  at  the  -center  of  the  hole,  rather 
it  is  the  pressure  existing  at  a  point  halfway  along  the  hole  radius  upstream  of 
its*  center.  Therefore,  the  correct  angle- between  the  centers  of  the  two  holes 


is  2  a 


+  180  d/?D.  where  a...  is  the  angular  location  of  the  point 


on  the  cylinder  where  the  pressure  is  just  equal  to  the  approach  flow  static 
pressure  and'd  and  D  are  the  diameters  of  the  hole  (0.020-irich)  and  of  the  cyl¬ 
inder  (0.083-inch),  respectively. 


+Static  pressure  obtained  using  the  static  pressure  holes  of  a  Pitot  tube  are 
usually  top  low  due  to  the  turbulence  velocities  normal  to  the  tube. 


Figure  6.  Schematic  of  Cylindrical  Pressure  Probe 


Flow 


Pressure  taps 


Figure  7.  Probe  Cross  Section  Showing  Pressure  Taps  and 
the  Angle  a 


25 


The  potential  £l<arsoluticn  for  the  flow  past  a  cirtular  cylinder  ih  cross 
flow  is  givenrbyAu  »..20SB  sin  a.  Apply  Bernoulli’s  theorem,"  the.  pressure  is 
given  by  p  *  p*,*  p(U„2  -u2)/2.  The  pressure  p  will  equal  the  uidisturbed 
static  pressure  whenu  =■  30?.  Thus,,  the  simplified  approachof  potential  theory 
gives  «  ^  *  -30®.  From  calibrations  in  steady  flow,  Glaser  [55]  fotnd  that 
“exit  “  measured ■acr.t  =  34.2°  at  Rep  =  8;5  x  103  arid  34.3° 

at  Rep  *  1.7  x  10*’.  Using  Glaser’s1  result  in  Thom's  equation  yields  = 

34.8  +  (90)  (0.020)?  (.083)  =  41.7°  for  the  present  probegeometry. 

The  slope  of  the  potential  flew  pressure  distribution  (£„  -  p)/pu2/2  is 
0.05  per  degree  near  a  pressure  coefficient  of  zero.  Therefore,  it  is. necessary 
to  determine  . 2a  to  within  0.4  degrees  in  order  to  prevent  a. systematic  error  of 
one -percent  in  1/2  pu2.  Since  the  slope  of  the  static  pressure  is  nearly  linear 


near  a^(  f  a  good  first  approximation  to  the  local  static  pressure  (or  flow 
angle):  can  be  obtained  by  taking  the  average  pressure  of  the  two  holes.  In 
reality  an  elaborate  calibration  was  required.. 

The  pressure  distribution  along  the  forward  surface  of  a  long  Pitot  cylinder 
varies  only  slightly  when  the  axis  of  the  cylinder  is  no  longer  set  perpendicular 


to  the  flow.  This  result  was  anticipated  by  Glaser  from  potential  flow  theory-- 
that  the  form  of  the  pressure  distribution  about  an  infinite  cylinder  is  indepen¬ 
dent  of  the  angle  between  the  free  flow  and  the  axis  of  the  cylinder.  Therefore , 
the  flow  can  be  broken  up  into  components  along  and  normal  to  the  cylinder,  the 


former  having  little  effect  on  the  pressure  distribution  since  the  velocity 
potential  has  the  property  c>f  being  additive  for  superimposed  flows.  In  the 
present  experimental  situation,  the  probe  is  always  aligned  normal  to  the  local 
total  velocity  while  the  crossflow  arises  because  of  the  swirl. 


Lavan  [16]  used  a^  =  40°  and  calibrated  P^W/P^g 


as  a  function  of 


Mach  number.  A  strong  Mach  number  dependence  on  the  pressure  distribution  was 
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observed.  Bruuh~[56]  alsd^calibrated.a  short  .Pitot  cylihder'in'a  free  jet 
(th.-mininrize,  wall  effects) .  Ihe.  observed  “jef  -f^pioxiaiately  45°)  was  constant 
for  velocities  under  100  m/sec.  Above  100-in/sec,  are£  has  increased  by  about 

finally,  Hinze  : [SI]  cautions  that  the  diameter  of-fbe  cylinder  imist  be 
•sniallj  at  least  hot  much  larger  than  the  microscale  of  the  ambient  turbulence. 

In  order  to.  measure  the  flow  angles  directly,  the  probe  was  firet  aligned 
such  that  >in  pure  ajtial  flow,  a  flew,  angle  reading  of  $  »  0°  is  obtained.  Values 
of  were  obtained  by  rotating  the  pressure  probe  mtil  the  indicated  static 
pressure  agreed  with  that  measured  using  a  static  wall  tap  located  in  the  inner 
tube  at  the  same  .axial-  (x/D^)  location.  Sample  calibrations  performed  in 
straight  flow  at  three  values  of  x/D^  are  given  in  Fig.  8.  variations  of 
±1/2°  are  observed  at  a  fixed  radial  location.  In  actual  practice,  each 
profile  was  plotted  on  a  larger  scale  and  a  smooth  curve  drawn  through  the  data, 
“ref  P10^1®5  at  intermediate  values  of  x/D^  were  obtained  from  interpolation 
of  the  calibration  profiles. 

Laufer  [57]  has  suggested  that  in  the  vicinity  of  the  wall,  velocity  correc¬ 
tion  should  be  made  because  of  the  large  local  turbulent  intensities  .  If  U  is 


the  total  velocity,  he  suggests  using 


U 


corr 


me  as 


K  urr  +  +  wtT 

/\1 - ^ - 


(16) 


For  normal  turbulence  intensities,  the:  correction  is -of  the  order  of  five 
percent  or  less.  Corrections  were  not  made  to  the  present  data  because  turbu¬ 
lence  intensities  were  not  measured. 

EXPERIMENTAL  PROCEDURE 

One  of  the  primary  objectives  of  this  study  was  to  obtain  velocity  profiles 
at  various  axial  stations  (D^  locations)  using  different  swirl  rates  (inlet  vane 
settings).  Velocity  surveys  were  obtained  at  stations  3,  7,  11,  15,  19,  24,  and 
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28  (x/D^  values  of  1.7,  4.2,  7.0,  10.3,  14.8,  22.2 ,  and  32.7)  for  nominal  Inlet 
vane  settings  of  0%  10°,  30°,  45%  and  60°.  The  bulk  of  the  data  involved - 
local  measurements  of  the  flew  angle,,  the  total  pressure,  the  static  pressure,: 
and  the  static  tenperature  of  the  flow.  All  of  these  quantities  are  heeded  to' 
obtain  the  local  velocity  conponents  at  a  single  point.  The  first  readings 
were  always  taken  at  a  radial  distance  of  0.010-inch  from  the  inner  tube  wall. 
In  all,  readings  were  taken  at  thirty  different  radial  locations  which  spanned 
the  annular  gap.  At  any  given  r-iocation,  the  flow  angle  was  determined,  first 
by  using  the  manometer  setup  as  a  U-tiibe  manometer  and  rotating  the  probe  until 
a. null  reading  was  obtained.  The  flow  angle  was  read,  to  the  nearest  five 
minutes ,  directly  -from  the  scale  on  the  base  of  the  traversing  mechanism.  Once 
the  flow  angle  had  been  determined,  the  probe  could  be  tinned  to  the  proper 
angles  for  measuring  the  total  and  static  pressure.  All  velocity  measurements 
were  made  at  sufficiently  low  velocities  (less  than  200  ft/sec)  that 
compressibility  effects  coiild  be  neglected. 

Tiie  axial  and  tangential  velocity  conponents  were  calculated  from  the  total 
velocity  vector  using  the  relations 


u  =  V  cos  <t> 

(17) 

and 

w  =  V  sin  <J) 

(18) 

where 

=  flow  angle 

u  =  axial  velocity 

< 

w  =  tangential  velocity 

A  sketch  showing  the  relative  position  of  the  total  velocity  vector  and 
the  flow  angle  <p  is  shown  in  Fig.  9. 
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Fig.  9 

Sketch  Showing  Flew  Angle  <}> 


The  mass  flow  rate  .v:as  determined  by  numerically  integrating  the  axial 


velocity  profiles.  Recall  that  with  swirl,  the  orifice  was  never  calibrated 
directly  to  obtain  the  mass  flow  rate.  The  mass  flow  rate  is  given  by  the 
following  integral 


m  =  2JI 


purdr 


(19) 


The  numerical  integration  was  carried  out  in  the  following  manner.  A  least 
squares,  second-degree,  polynomial  was  fit  through  the  first  four  data  points, 
with  the  first  point  how  being  taken  as  the  zero  velocity  point  at  the  inner 
wall.  The  least  squares  fit  was  of  the  form  ur  =  Cj  +  C2r  +  C3rz.  This 
polynomial  was  then  integrated  between  the  wall  and  the  third  data  point,  th;is 
giving  the  mass  flow  rate  through  that  interval  as 


Ai  =  2IIp{Ci[r(3)  -  ri]  +  C2/2[r(3)2  -  r,2]  +  C3/3[r(3)3  -  r^]}  (20) 

The  local  density  used  was  the  average  of  the  densities  at  the  points  over 
which  the  integration  was  performed.  Away  from  the  inner  wall,  a  set  of  five 
points  was  used  for  the  least  squares  fit  and  the  integration  was  carried  out 
between  the  second  and  third  point  of  the  five -point  set.  The  first  point  of 
each  set  was  the  same  as  the  second  point  of  the  previous  set. ,  At  the  outer 
wall,  the  fit  was  again  through  four  points  and  the  integration  was  carried  out 
over  two  intevals. 


The  radial  velocity  gradients  were  alsocalculatedat  each  point-  Essen¬ 
tially,  the  sane  technique  was  used  to  differentiate  the  velocity  profiles  as 
was  used  to  integrate  them.  Instead  of  integrating  a  fitted  least  squares  poly¬ 
nomial,  a  differentiation  was  required.  For  exaaple,  the  least  squares  ^poly¬ 
nomial  fitted  to  the  axial  data. was  of  the  fbim  u-'=  Gi  +  C2  *  y  +  C3  *  y2, 
where  y  is-  equal  to  the  radial  distance  measured  from  inner  wall  (r  -  r^). 

The  velocity  gradient  was  evaluated  at  the  center  point  of  each  five-point 
fit.  At -each  wall,  a  four-point  fit  was  used  and  the  gradient  was  evaluated 
'at  the  wall  and  the  next -two  points. 

DISCUSSION  OF  RESULTS 

Axial  Velocity-  Profiles.  The  extreme  axial  velocity  profiles,  presented  in 
Figs.  10  and  11,  were  plotted -at  the  initial  and  final  axial  location  with 
inlet  vane  angle  as  the  parameter.  A  somewhat  unexpected  profile  exists  at 
station  3  (Fig.  10).  For  zero  swirl,  the  axial  velocity  profiles  have  a  nearly 
felly-developed  characteristic  shape  although  located  quite  close  to  the 
entrance  (x/I^  =  1.7)  where  centrally  flat  profiles  were  expected.  It  appears 
that  the  friction  drag  of  the  entrance  plates,  coupled  with  the  profile  drag 
of  the  turbulence  trips,  has  strongly  retarded  the  flow  near  both  walls,  there¬ 
by  producing  the  more  rounded  velocity  profiles  characteristic  of  fully  dev¬ 
eloped  turbulent  annular  flows.  Further  downstream,  the  calculated  average 
axial  velocities  of  the  $  =  0°  -*■  45°  data  are  constant  to  within  +3 .percent, 
while  the  <jr=  60°  profiles  seem  to  have  changed  considerably  in  comparison  with 
the  other,  profiles-. 

The  axial  flow  field  is  most  sensitive  to  inlet  vane  angle  $  greater  than  45° 
For  J  -  60°  in  the  inlet  region,  the  flow  is  approaching  a  reversed  flow  condi¬ 
tion,  (axial  velocity  near  the  inner  wall  moving  in  the  negative  x-direction) , 
which  explains  the  skewed  nature  of  the  60°  profile.  Since  the  inner  flow  is 


Velocity  Profiles.  ‘Pressure  Probe,  Station  3 


Velocity  Profiles .  Pressure  Probe ,  Station  28 


TCte'nie^i  -the'.duter.  fldw.must.be :  accelerated, to. a  larger  than.^ormal -velocity 
to  cofipehsatei  F^ow  reversal  -did-.occur  at.  an  inlet  Vane  setting.of  approximately' 
$-=  65/-  jthe  calculations  of  .Bossel  .  [4],  •  The  $  -=  $0°  profiles 

developed  in  a  regular.manner,  becoming  flatter  and  fuller. -as  one  moves  down¬ 
stream..  .This-  data:  ejdriibited  mpre.  msteadiness  and  larger  than. normal  scatter-  • 
resulted. 

*  f 

For  straight  flow,  the  maximun- axial  velocity  at  station -3  is  located  at  a 
dimensibhless  radial  position  -(r  -  T^)/(r0  -  r.)  of  about  .0.46.  Downstream  of 
station,  3  the  velocity  maximuiufirst  moves  closer  to  the  inner  wall  and  then  out- 
ward  to:;a  final; yalue  of  0; 45  at  station  28.  .For  <j>  *  10°,  the  maximum  location 
begins  at  0:37  at  station  3,  moves  inward-  and  then  outward  until -at  station  28 
it  occurs  at  Arm/ArQ  «  0.50.  The  same  trends  are  noted  for  J  =  30°  and  "ij>  =  45°. 
The  values  reported  are  .only  approximate  since  they  were  obtained  by  visual  " 
inspection  of  the  graphs,  while  the  irregularities  in  the  location. of  the 
-maximum  velocity  are  clearly  observed.  For  b  3  60°,  the  maximun  velocity  occurs 
much  closer  to  the  outer  wall.  At  station  3,  the  maximun  occurs  at  a  dimension¬ 
less  radius  of  approximately  0.8  and  does  not  vary  from  this  location  as  the  flow 
proceeds  downstream. 

Tangential  Velocity  Profiles.  The  tangential  velocity  profiles  obtained 
using  the  pressure  probe  are  plotted  for  each  inlet  vane  setting  with  the  station, 
i.e.,  dimensionless  axial  location,  x/D^,  as  -the  parameter.  The  axial  decay  of 
tiie  .tangential  velocity  component  can  be  readily  observed  by  plotting  the  profiles 
in  -.this  manner.  The  profiles  will  be  discussed  starting  with  those  for  ?  =  10°. 
Fig.  12. 

One  obvious  characteristic  of  the  profiles  is  the  existence  of  two  relative 
maximum. points  for  over  half  of  the  profiles  at  <j>  =  10°.  The  swirl  is  also 
relatively  weak  with  the  absolute  maximum  tangential  velocity  at  station  3  being 
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.only  15%  "of  the  corresponding  maxiiua  axial  velocity.  ;The-  absolute  soxinua 
of  e^h- profile-  ;6ccip,s  close  to  the  outer  wall  ,  i.e.,,  ata-dimensianless  radial 
position. of, Of 9.  The  smaller  relative  maximun  points,  when  they  occur,  are 
hear  co  the  inner  wall  (around  0.3);  further  downstream,  the  existence  of  the 

two  relative  maximum  points  becomes  less  apparent,  until  at  station  24,ohiy 

*v  '  - 

one  maximum  point  is  £om3.  ■  . 

The  flow  adjacent  co  the  inner  boundary  layer  is  similar  to  a  solid  body 
rotation  (forced  vortex  type)  as  is-  evident  at  stations  24  and  28.  Closer  to 
the  inlet  the  forced  vortex  flew  gradually  reverts  to  a  free  vortex  form  which 
extends  outward  toward  the  outer  wall.  Adjacent  to  the.  outer  wall  the  flow 
reverts,  back  .to  the  forced  vortex  form,  and  produces  the  double  maximun. 

At  ^  =  309,  the  maxut>un  tangential  velocity  at  station  3  is  nominally  50% 
of  the  maximun  axial  velocity.  The  existence  of  two  relative  maximum  points  is 
noted  only  at  stations  3,  7,  11,  and  15.  At  station  3  the  absolute. maximum 
point  is  at  a  Arm/Aro  of  0.13.  The  maximum  points  at  stations  7,  11,  and  15 
are  again  close  to  the  outer  wall.  One  now  observes  that  at  station  3  the 
majority  of  the  flow  is  of  a  free  vortex  type.  At  4>  =  30°  the  vortex  generator 
is  operating  as  planned.  It  is  of  interest  to  note  the  rapid  change  from 
essentially  a  free  vortex  type  in  the  outer  2/3  of  the  annulus  at  station  15, 
to  a  completely  forced  vortex  type  at  station  19.  The  forced  vortex  pattern 
is  even  more  evident  at  stations  24  and  28. 

At  $  -  45°,  Fig  14.,  the  maximum  tangential  velocity  is  now  more  than  twice 
that  observed  at  $  =  30°.  The  maximum  tangential  velocity  at  station  3  is  103% 
of  the  corresponding  maximum  axial  velocity.  The  profile  at  station  3  is  the 
only  one  that  exhibits  two  relative  maximum  points.  All  of  the  maximum  points 
occur  within  the  inner  half  of  the  annulus  except  at  station  28  where  the 
maximum  occurs,  at  a  dimensionless  position  of  about  0.62.  All  of  the  profiles 


exhibit  a.free  vortex  type  of  flow  over  a;majpr  portion  of  the  annular  flow 
regibn:  and'the  See,  vortex  nature  is  still  preserved  at  station  28.  "Die  free 
vortex  produced  by  the  inlet  guide  vanes  isbecoedng  stronger,  i.e.,  a  larger 
total  angular'  noBEntum  is  initially  produced  and  is  making  its  presence  felt 
over  the  entire  length  of  the  test  section.  The  forced  vortex  type  of  flow  is 
still  evident  close  to  the  inner  wall  and  itmoves  out  until  at  station  28  the 
flow  in  the  inner  half  of  the  annulus  has  a  forced  vortex  character. 

At  ^  =  60°,  Fig.  15,  the  aiaxinuh  tangential  velocity  at  station  3  is  new 
approximately  160%  of  the  corresponding  maximum  axial,  velocity.  No  double- 
maximtm  profiles  occur.  All  of  the  profiles  new  exhibit  a  maximum  point  at  a 
dimensionless'-positicu  i(within  10%)  of  6.5,  with  the -maximum  at  station  28 
being  closer  to  the  inner  wall  than  the -..aximum  at  station  3.  The  annular  flew 
region  is  thus  approximately  divided  in  half  with  the  inner  half  exhibiting 
forced- vbttex  flow  and  the  outer- half  exhibiting  free  vortex  flow. 

To  sunmarize,  the  tangential  velocity  curves  for  each  inlet  vane  setting 
reveal  the  following  regions: 

1)  Inner  wall  boundary,  layer  produced  by  the, no-slip  boundary  condition  at 
the  inner,  impermeable  wall,  i.e.,  at  r  =  r- ,  u  =  v  =  w  =  0,  dw/dr  >  0. 

2)  Forced  vortex  zone.  The  region  above  the  slight  discontinmty  in  shape 
of  the  swirl  velocity  profile.  The  tangential  velocity  increases  with 
increasing  radius  according  to  the  approximate  relation. w  =  const  r*11 , 

0  <  m  £  1  and  the  constant  is  a  function  of  axial  position- x. 

3)  Transition  zone.  The  region  brackets  the  location  of  maximum  tangential 
velocity. 

4)  Free  vortex  zone.  The  tangential  velocity  is  described  by  approximate 
relation  w  =  const  rm,  -1  <_  m  £  0. 

5)  Outer  wall  boundary  layer  occurs  in  the  region  adjacent  to  the  outer  wall 
The  swirl  velocity  gradient  is  negative  and  at  r  s  r  ,  u  =  v  a  w  =  0. 
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The  static  pressure  distributions  along  the  inner  wall  are  presented  in 
Fig..  16.  Since  atmospheric.. pressure  (Patja)  exists  at  the  inlet,  all  of  the 
measured  pressures  are  siiatmospHeric.  The  pressure  differences-,  at  . x/I^  =  0 
represent  the  entrance  pressure  losses  and  are  due  primarily  to  vortex  gener¬ 
ator  friction;  The  flow  cross  sections  within  the  vortex  generator  become 
quite  small  at  =  60°.  The  airfoil  sections  nearly  overlap.  The  average 
velocity  within  the  vortex  generator  increases  sharply. 

The  slope  of  the  profiles  is  given  by 


3(P0._  -  P)  .w 

-ja*S _ _  -n  12_ 

3  x/T^  n  3x 


(21) 


Hie  normal  behavior  of  a  decreasing  static  pressure  with  axial  distance 
is  observed  for  ^  =  0°,  10°,  and  30°.  An  extensive  region  of  nearly  uniform 
pressure  is  observed  in  the  region  6  <  x/D^  <  22  for  $  =  45°.  The  internal 
wall  pressure  gradient  is  positive  over  the  entire  test  length  for  f  =  60°.  There¬ 
fore,  tire  normal  physical  feeling  that  there  can  be  no  extensive  regions  of  flow 
against  a  pressure  gradient  is  clearly  not  true  in  rotating  flow  systems. 

Further  Data  Reduction  and  Tests.  Mare  extensive  discussion  of  the 
$  =  45°  data  will  be  presented,  here  because  this  data  represents  the  results  of 
an  intermediate  swirl.  The  axial  velocity,  flow  angle,  tangential  velocity, 
and  angular  momentum  profiles  are  given  in  Figs.  17,  18,  19,  and  20,  respec¬ 
tively. 

Local  Skin-Friction  Coefficient.  Local  axial  skin- friction  coefficients 
were  computed  for  both  inner  and  outer  walls.  The  tniversal  velocity- distribution 
equation  for  turbulent  boundary  layer  flew  has  the  form 

u+  *  A  •  log  y+  +  B  (22) 


Static  Pressure  Distribution  Along  Inner  Wall 
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where 

uv:  r  u/u  (23) 

+ 

y  =  y  •  u  /v  (24) 

Ihe  values  of  tie  law  of  the  wall  constants  A  and  B  vary  somewhat  as  is 
.  observed  in  the  following  table.  .  5 


Table  I.  Law  of  Wall  Constants  1 


Source 

A 

B 

Nikuradse 

s.s: 

5.8 

Clauser 

5.6 

4.4 

Coles 

5.62 

5.0 

Smith  and  Walker 

5.0 

7.15. 

Patel 

5.5 

5.45 

NPL  Staff 

4.9 

5.9 

B  assumes  different  values  for  roughness  as  summarized  in,  [58],.  Brighton  and 
Jones  ;[59]  have  carried  out  careful  experiments  in  annuli  with  several  radius 
ratios.  They  find  the  values  of  ii+  generally  lie  close  to  the  distribution  of 
Nikuradse  [58]  for  the  outer  profile,  as  do  the  results  of  Knudsen  and  Katz  [60], 
The  inner  profile  results  for  r-/rQ  =  0.56  agree  better  with  the  Clauser  constants. 

Similarity  considerations  by  Eskinazi  and  Yen  [11]  describe  deviations  from 
the  universal  law  for  smooth  walls  due  to  the  radius  of  curvature  of  the  channel. 
For  convex  walls  (inner  wall)  the  velocities  are  larger  than  predicted  by 
equation  22,  whereas  for  concave  (outer  walls)  the  velocities  are  less  than  pre¬ 
dicted  by  the  straight-wall  universal  law. 

The  values  of  A  and  B  used  herein  were  those  of  Nikuradse.  Using  these 
constants,  the  equation  22  represents  the  turbulent  core  region,  i.e.,  y+  >  30. 
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'Makinguse  of  the  fact  that 

u*M  e  (0^2) 1/2 

where 

u*  =  /r  /p  ;  axial  shear  stress  velocity,  ft/sec 
wi 

ui  =  axial  velocity  at  outer-edge  of  inner  axial  boundary 
layer,  ft/sec.  '(In  this  case,,  it  Is  also  equal  to  - 

,  .the  maximum  axial  velocity.) 

Cf  *  Tw  /^Pui2  »  local  axial  skin-friction  coefficient 
tx  wi  at  inner  wall 

Equation  22  can.be  rewritten 

u/ui  =  (cfx/2)1/i  [2.5,  In  (yui/v).  +  2.S  In  (Cf/2)1/2  +  5.5]  (27) 

Thus,  tliis  equation  gives  a  curve  of  u/ui  versus  y*Ui/v  for  .each  value 
cf  the  parameter  C£x  (see  Fig.  21)  .  As  suggested  by  Clauser  [61],  values  of 
C£x  now  be  estimated  by  plotting  the  experimental  values*  of  u/ui  versus 
y  u\/\)  and  selecting  the  Cfx  value  that  best  fits  the  data.  Equation  27  was 
plotted  on  semi-log  paper  so  that  curves  of  constant  Cfx  would  plot  as  straight 
lines. 

The  values  of  Ui  needed  to  determine  C£x  were  determined  from  the  appro¬ 
priate  axial  velocity  profiles.  For  the  axial  velocity  profiles,  the  value  of 
Ui  was  ‘-the  same  as  the  maximum  axial  velocity. 

Pierce  [62]  has  given  a  good. discussion  of  a  Clauser  plot, as  applied  to  a 
skewed  boundary  layer  flew.  At  least  four  variations  of  the  conventional  Clauser 
plot  ate  possible  for  a  three-dimensional  turbulent  flow.  He  compares  experi¬ 
mental  results  of  pressure  driven  skewed  turbulent  boundary  layer  flows  using 
calibrated  fteston  tubes  [63],  Clauser  plots,  claw  Preston  tubes,  a  direct 
force  balance  system,  and  a  wall  heat  meter.  Pierce  concluded  from  accunulated 
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Universal  Velocity  Distribution— Inner  Profile 


so 


data  ;that  the  Clauser  chart  described  above  should  predict  local  wall  shear 
values  with  as  much  accuracy  as  in  the  two-dimensional  case  provided  a  reason- 
’■  able  nuntoerofmeasured  velocity  points  are  recorded  in  the  wall  similarity 
region. 

The  experimental  data  plotted5 in. Fig.  22  represents  some  of  the  best  fits 
arid  the  two  worst  fits  of  the  <jT  =  -459  data  to  lines  of  constant  C^.  For  each 
flow  angle,,  stations  3,  11,  10,  and  28  generally  represent  the  worst  to  the 
best  fits,  respectively; 

Hie  actual  plot  used  to,  deteimine  the  values  of  C£X  has  an  enlarged  linear 
scale  of  10  squares  per  inch-,  with  one  inch  representing  an  increment  of  0.02 
in  u/ui,  and  an.  enlarged  log  scale  with  one  cycle  covering  approximately  10 
indies.  Lines  of  constant  were  drawn  in  at  increments  of  0.0002.  A  linear 
inteipolation  was  used  between  these  lines  to  determine  the  values  of  the  constant 
lines  that  were  drawn  through  the  data  points.  For  data  points  that  did 
not  fail  along  lines  of  constant  Cfx,  a  straight  line  was  drawn  through  the  data 
points  and  the  value  of  was  evaluated  at  the  point  where  a  line  of  constant 
y*Ui/\)  =  5,000  intersected  the  line  through  the  data  points.  The  method  is 
illustrated  in  Fig.  21  for  station  3,  $  =  45°. 

Law-of-the-Wall  Variables.  Once  the  axial  skin- friction  coefficients  had  been 

+  + 

determined,  the  law-of-the-wall  variable  for  the  inner  profile,  u  and  y  ,  could 
be  calculated.  Making  use  of  equation  25,  the  values  of  u+  and  y+  can  be  easily 
calculated  from  their  definitions. 

Hie  logarithmic  velocity-distribution  equation  for  turbulent  flow  is  attri¬ 
buted  to  Prandtl,  who  used  his  mixing  length  theory  as  a  basis.  As  discussed  in 
the  previous  section,  the  constants  used  in  the  universal  velocity  law  were  those 
of  Nikuradse.  Using  his  values  for  the  constants,  the  equation  becomes 

u+  =  2.5‘ln  y+  +  5.5 


(28) 
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Values  of.  u+  versus  y+  were  then  plotted'  oh  semi-log  paper  (see  Fig.  22)  and 
conpared  with  the  assumed,  relation  given-by  equation  28. 

*«r  v  '  ^ 

Dimensional  values  of  the.  axial  and  tangential  wall  shear  values  are 
plotted  in  Figs.  23  and' 24.  Inner  wall  axial  shear  values  obtained.. separately 
from  a  calibrated,  Preston  tube  and  the  Pitot  cylinder  probe, -Clauser  chart  are  in 
agreement.  Outer  wall  tangential  shear-  -ess  values  are  approximately  twice 
as  great  as  the  inner  wall  values,  reflex  ig  the  large  tangential  velocity 
gradients  near  the  outer  wall  shown  in  Fig.  19.  The  inner  wall  shear  data 
exhibit  a  wavy  behavior  in  the  range  x/D^  <  10  which  reflect  the  pressure 
distributions  plotted  in  Fig.  16.  Axial  pressure  distributions  are  given  for 
(r  *•  fp/Cfg  ■  r-)  »  0  (inner  wall),  and  0,2,  0.4,  0.6,  0.8,  and  1.6  (outer  wall). 
For  x/D^  >  10,  the  inner  wall  pressure;  gradient  is  negative  but  small.  At 
other  radial  locations,  the  local  pressure  gradients  are  positive. 

The  entrance  pressure  losses  are  largest  near  the  outer  wall.  Here,  the 
pressure  gradient  inthe  axial  direction  is  due,  on  one  hand,  to  the  effect  of 
accelerating,  and  overaccelerating,,  the  flow  (which  initially  has  no  axial  com¬ 
ponent)  to  the  velocity  it  attains  at  station  x.  On  the  other  hand,  ultimate 
radial  pressure  equilibrium  and  the  axial  decay  of  angular  momentum  lead  to 
complex  pressure  patterns. 

Local  axial  shear  coefficient  profiles  for  <F  =  45°  are  given  in  Fig.  26  at 
five  axial  stations .  These  values  were  obtained  by  numerical  integration  of 
equation  9',  omitting  the  turbulence  intensity  change  term  (term  5).  The 
radial  location, for  zero  axial  shear  moves  outward  from  Ar/ArQ  =  0.34  at  station 
7  to  Ar/ArQ  -  0.61  at  station  24.  The  axial  velocity  profiles  attain  their 
maximum  values  (ui)  at  Ar/ArQ  =  0.52  at  station  7  and  0.58  at  station  24.  This 
discrepancy  will  prove  to  be  serious  when  eddy  diffusivities  are  computed  later. 
The  outer  portions  of  the  profile  at  station  7  indicate  a  thick  constant 
shear  region. 


Pressure  Distribution  at  Several  Radial  Locations 
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Contribution  of  the  individual  terms,  in  equation  9',  which  sum  to  the 
local  shear,  are  given  in  Fig.  27  for  ^  =  45°,  at  station  7.  Term  1  (circles) 
represents  the  measured  inner  wall  shear  modified  by  the  radius  ratio  r^/r.  Term 
2  (triangles)  involves  axial  derivatives  in  the  integral  of  the  axial  velocity. 
Tiiis  second  term  is  zero  at  each  wall  and  nearly  symmetrical  because  of  axial 
velocity  near-syumetry.  At  the  outer  wall,  the  integral  represents  the  average 
axial  velocity  which  dees  not  change  significantly  with  x.  As  the  flew  develops, 
the  velocities  decrease  such  that  term  2  is  positive.  Term  3  (inverted  triangles) 
result  from  changes  in  kinetic  energy.  These  points  are  all  negative  and  nearly 
symmetrical  about  the  radial  location  for  maximum  axial  ■velocity.  The  fourth 
term  (squares)  represents  the  net  pressure  force  term.  This  term  is  very  small 
near  the  inner  radius  but  grows  due  to  the  increased. swirl  near  the  cuter  wall. 

The  total  shear,  plotted  as  the  x’s,  is  the  sum  of  terms  1-4. 

Similar  conparisons  are  available  in  Fig..  28.  Here  the  conposition  of  the 
axial  shear  at  Ar/ArQ  =  0.2  is  given  at  5  axial  stations.  This  radial  location 
was-.selected  because  the  swirl  velocities  are  normally  near  their  maximum  values 
nearby.  This  figure  does  not  answer  the  question  of  exactly  how  the  inner  wall 
shear  obtains  its'  value.  It  is  simply  a  measured  value  and  the  integral  equation 
9'  demonstrates  why  the  shear  changes  radially.  The  effect  of  swirl  enters  only 
indirectly  through  the  pressure  and  the  axial  profile  modifications.  At 
Ar/Ar^  =  0.2  the  acceleration  and  kinetic  energy  terms  tend  to  cancel  each  other 
put  while  the  pressure  term  remains  small  such  that  .the  local  shear  is  nearly 
equal  to  the  inner  wall  value. 

Conparable  local  tangential  shear  curves,  resulting  from  the  numerical 
integration  and  differentiation  of  equation  11',  are  given  in  Figs.  29  and  30 
fer  f  =  45°.  Near  the  inlet  xr0  “  0  at  very  small  Ar/ArQ  values  (less  than  0.1). 
At  station  24,  x^  ~  0  at  Ar/Ar0  “  0  •  4^1  •  Near  the  inner  wall,  xrQ  is  propor¬ 
tional  to  Ar/ArQ  whereas  near  the  outer  wall  xrQ  is  constant. 
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The  magnitudes  of.  the  individual,  teims  in  equation  11*  .are  plotted  in  Fig.  30 
for  ^;=  45°',  station  7.  -  The  circles  represent  term  1,  tie  torque -per  unit  length 
on  the  inner  wall  modified1  by  a  radius  ratio  squared.  Terms  2  and  3  are  the 
change  in  axial  flujcof  ..angular  momentum,  and  the  effect  :of  acceleratidri  on  the 
axial  flux -.of  angular  momentum.  These  terms  are  of  opposite,  sign  arid  tend  to 
cancel.  Symmetry  about  the  position  of  maximum  axial  velocity  is  weak  because 
the  swirl  velocity  is  quite  unsynmetric.  Tne  tangential  shear  changes  most 
rapidly  in  the  radial-  direction  near  the  zero  shear  location. 

The  dominant  term  is  the  change  in  axial,  flux  of  angular  momentum  between 
the  inner  wall  and  location  r.  The  shear  "bulge"  located, in  the  range  0.2  < 
Ar/Ar0<  0.8  appears  to  be  an  effect  of  the  entrance  condition  which  has  disap¬ 
peared  between  station  IS  and  19.  The  tangential  shear  profiles  at  stations  19 
and. 24  are  quite  similar  in  shape,  indicating  that  the  entrance  effects  are 
probably  finished  aid  the  final  process  of  swirl  decay  is  underway. 

The  tangential  shear  near  the  outer  wall  is  approximately  twice  as  great 
as  the  inner  wall  sheaf.  This  reflects  the  existence  of  larger  swirl  velocities 
near  the  outer  wall  and,  presumably,  the  increased  turbulence. 

Axial  variations  of  the  shear  conponent  terms  are  plotted  in  Fig.  31  for 
location  Ar/ArQ  =  0.2,  ^  =  45°.  In  the  initial  swirl  decay  region,  x/D^  <  12, 
the  local  shear  is  determined  almost  entirely  by  the  change  in  the  axial  flux  of 
angular  momentum.  As  the  swirl  decays  and  becomes  progressively  weaker,  a  type 
of  fully  developed  situation  arises  where  the  two  integral  teims  in  equation  11’ 
approximately  cancel  while  both  tend  to  zero,  such  that  the  local  shear  is 
nearly  equal  to  the  inner  wall  shear. 

EDDY  DIFFUSIVITIES 

Axial  Diffusivities.  Axial  eddy  diffusivities  were  computed  from  the  relation 
*rx  1? 
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The  a^ial  diffusivities  are  plotted  in'  Figs.  32-36  as  the  ratio  of  turbulent 


v\ 


to  molecular  diffusivities  plus  uiity  ’(erxA>  +  1)  vs  the  radial  station  Ar/ArQ. 
The  ordinate  is  the  ratio  of  the  total  diffusivity^to  the  molecular  diffu- 

_  \  *  •  i 

sivity  for  air  at  the  same  tenperature.  The  axial  diffusivity  data  are -lowest 
near  the  walls  in  all  cases.,  Diffiisivi 


sivities  for  Ar/ArQ  >0.5  usually  exceed 
o 

same  radial  distance  from  the  inner  wall  as.Ar/Ar^  =0.9  is  from  the  outer  wall. 


corresponding  values  below  0.5.  For  example,  a  Ar/Ar0  of  0.1  is  located  the 


At  station  7, 'that  data  for  zero  swirl  are  represented  by  the  circles.  One 


point  neaij'  Ar/Ar  =  0.4.  is  missing  because  the  shear  stress  was  opposite 


in 


sigh  to  the  velocity  gradient.  The  observed  diffusivity  ratio  exceeds  600 


>jr  4>  =  0. 


roar  the  inner  wall  for  $  =  60°  as  conpared  with  a  voice  of  120  fo| 

A, great  many  60°'  swirl  joints  are  missing,  near  .the  outer  wall.  The  apparently 
"wild" -points  in  the  central  areas  of  the  figurJ  are  misleading.  All  of  the 
data- is  systematic.  The  apparent  scatter  is  associated  with  the  onset  of  zero 
sheaf  arid  zero  gradients.  The  total  axial  diffusivity  ratio,  as  confuted  from 
measured  shear  stress  and  measured  velocity  gradient,  goes  to  infinity— plus 
or  minus— in  a  variety  of  ways.  Mast  of  the  effects  observed  in  the  grange 


0.3  <  Ar/ArQ  <,  0.7  are  due  to  this  effect. 


I 


*  I 

The  axial  velocity  profiles  are  nearly  fully  developed  by  station  24. 

Therefore,  comparisons  may  be  made  with  other  fully  developed  diffusivity 
results  and  predictions.  The  experiments  of  Jonsson  [64]  are  closely  appli¬ 
cable.  He  worked- with  a  concentric  annulus  with  a  ratjius  ratio  of  r^/rQ  =  0.56 
(conpared  with  the  present  value  of  0.4) .  His  diffusivity  data  apply  at 

x/Dh  ~  132  and  an  axial  Reynolds  number  of  Re„  =  115,000,  The  agreement  is 
\  ‘  “ 

fair  although jJonsson's  inner  data  lie  above,  and  the  outer  data  below,  the 


present  zero  swirl  dat4  at  station  24-. 


\ 


! 


min 


Figure  36.  Turbulent  Axial  Eddy  Diffusivity— X/E^a  22.2 


69 


Two  enperical  relations,  which  confute  a  constant  diffusivity  for  solid 
body  rotation  in  astationary  pipe,  .are  plotted  in  Fig.  36.  They  are  labelled 
Fejer  and  Javan  [16] 


(8. 32) (1.3x10s)0* 3  =  285 


(30) 


and  Krieth  and  Sonju  [37]* 


v 


5.44x10 


0.86 


V 


) 


=  (. 00544) (1. 3x10s) °* 86  =  135.9 


(31) 


Away  from  the  wall  boundary  layers,  most  of  the  present  data  at  station  24  fall 
within  these  values. 

Recently,.  Wassan,  Tien,  and' Wilke  [65]  pointed  out  that  most  of  the  pro¬ 
posed  (axial)  diffusivity  distributions  do  not  satisfy  the  theoretical  criterion 
of  Townsend,  which  requires  that  the  turbulent  contribution  to  the  Reynolds 
stress—  pu'v1  is  proportional  to  yn  near  the  wall  where  n  is  not  less  than 
three.  They  propose  the  relation 

=  4.16xl0~y3  -  15. 15x10' 6y+4 

v  1  -  4.16xl0’“y+3  +  15. 15x10' 6y+4  (32) 

which  applies  for  0  <_  y+  £  20. 


*  The  predictions  of  Krieth  and  Sonju  are  based  on  data  presented  in  an  M.S. 
Thesis  by  Musolf  [40].  Mussolf's  data  is  not  conplete  because  he  did  not 
measure  either  the  axial  or  the  tangential  velocity  profiles.  His  tangential 
velocities  were  calculated  from  measurements  of  the  local  flow  direction,  the 
total  mass  flow  through  the  pipe,  and  an  assumed  axial  velocity  distribution. 


*'-?  y  values  measured  in  the  present  study  are  found  in  Fig.  22  to  be 

20.  Bor  a  y+  of  20,  the  above  equation  yields  e TX/v  =  10. 

For  siqple  mixing  length  theories  inwhich  t  =  p£2(3u/3y)2,  £  =  ky 


formally,.*  =  0.4  and  s^/v  =  8  for  y+  =  20.  It  is  clear  from  these  values  that 
the  present  data:  were  not  taken  sufficiently  close  to  either  wall  to  distinguish 
any  viscous  layer  details. 


Tangential  Diffusivities.  The  tangential  eddy  diffusivities  were  conputed 
using  the  defining  relation 

he  =  pCv  +  \eHfj  -  f>  (3,) 

Again,  the  results  are  plotted  vs  Ar/ArQ  in  Figs  37-41  with  the  ordinate 
taking  the  form  of  the  ratio  of  total  tangential  diffusivity  to  molecular 
diffusivity  (erg/v  +  1) .  In  general  cr0  <  e^.  The  values  of  er0  are  quite  low 
near  the  inner  wall-particularly  at  the  large  swirl  rates.  There  is  a  definite 
decrease  in  erQ  with  $  near  the  inner  wall  at  station  7  and  at  most  of  the 
following  stations. 

The  enperical  predictions  of  Fejer  and  Lavan  [16]  and  Kreith  and  Sonju  [37] 
for  the  decay  of  solid  body  rotation  in  pipes  are  plotted  for  reference.  They 
represent  average  values  in  several  senses.  First,  they  apply  to  both  the  axial 
and  the  tangential  diffusivities.  Secondly,  they  do  not  vary  with  the  radius. 
The  predictions  agree  best  near  the  outer  regions  of  the  annulus  where  most  of 
the  annular  area  occurs. 

The  fully  developed  curved  channel  experimental  diffusivity  results  of 
Wattendorf  also  appear  in  Fig.  41.  Near  the  inner  wall,  the  straight  channel 


Figure  i*0.  Turbulent  Tangential  Eddy  Diffusivity— X/E^-  ilj.8 


\ 


'  .A  ^  \ 

\  .  \  1 

diffiisivities  are  approximately  one-jialf  the.  curved  wall  values.  Near  the 

outer  wall,  the  reverse' is  true.  The  locations  of  the  zeroes:.  3w/3r  -  0, 

*1,  * 

"  ,  **  ’  >}’*  » 

(3w/3r.-w/r)  --  and  (3w/3r.+  w/r)  0  closely*  match  the  .present  <(>  =  30°  data. 

In  fact,  the  similarity  in  behavior  in  all  respects  suggests  that  the 

Wattendorf  curved  channel  flow  and  the  present  flow  are  of  tlie  same  class.  ^ 

I  As  is  the  case  with  the  axial  momentum  transport ,  the  locations  of  zerr 

<shear  do  not  occur  at -the  same  radial  location.  -Infinities -,in;e0  are  appro? 

clied.  This  \is  the  same  type  of  "apparent!  scatter"  iLund  in  the  cxial  diffusivity 
*  | 

data.  In  general,  the\zero  tangential  shear  location. occurs  closer  to  the  inner 
wall  than  the  zero  axial  shear  because  of  the  subtraction  in  the  rate  of  strain 

.......  .  .  .  .  .  .  i 

term  (3w/3r  -  w/r).  The  confuted  diffusivities  do  increase,  reaching  maximum 
values  at  the  edge  of  the  tangential  boundary  layer  adjacent  to  the  outer  wall. 
Lavpn  and  Fejer  [16]  have  suggested  that  in  a  decaying  swirling  flow  in  a 

(  i  « 

station  pipe,  th^  eddy  diffusiyity  is  a  strong  (approximately  cubic)  function  of 
the  radial  coordinate.  This  trend  is  usually\ followed  by  the  present  data 
outside  of  the  point  of  zero  shear.  Finally,  the  diffusivity  values  in  the 

outer  swirl  boundary  layer  decrease  as  the  wall  is  approached.  • 

'  \  » 

Turning  to  the  la$t  staci'  ,  Fig.  ^1,  the  different  behavior  near  the  two 

i  ; 

walls  is  evident.  The  tangential  diffusivities  first  increase  with  Ar/Ar0. 
Outside  o|  this  "inner"  boundary  layer  -er0  decreases,  appjro  aching  a  value  of 
less  than  10. 

As  was  -anticipated  in  the  earlier  sections,  axial  diffusivities  near  the 

inner  wall  are  considerably  lower  than  corresponding  values  near  the  outer  wall.] 

-  1 

I 

This  is  a  res.ult  anticipated  in  the  Introduction  and  is  due,  presumably,  to  the 
different  turbulence  'characteristics  of  the  two  layers.  Also,  there  appears 
to  be  a  slight  increase  in  outer  wall  axial  and  tangential  diffusivities  with  1 
an  increase  in  swirl. 

Ait  alternate  form  of  presenting  the  dimensionless  diffusivities  is  revealed 


1 
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,\ 


in  Figs.  42-51.  The’  raticj  of  tuibulent  diffusivity  to  the.  product  oJ  wall 

friction  velocity  times'  ah  appropriate  =  scale  length .^.s  the  inverse  of  a  .Reynolds 

* 

number..  The  .flow  is  divided  into  two  parts  by  the  zero-shear  radius.  The  two 

\  . \ 

correlation  parameters  are  ; 


rx 


ui*CrT  ,.w»: 

rx  " 


or 


rx 


ri? 


”  uo*(r0 :  rt  J 

“i  rx  u 

i 


(35) 


where 


Uj*  =  inner  wall  friction  velocity  =  /x  '/'p 


lip  =  outer  wall  friction  velocity  =  /r^/d 

r_  n  = '(radial  location  for  zero  axial  shear' 
•'Trx=  0  '  I 


r^  *»  inner  wall  radius 
rQ  =  outer  wall  radius 


\ 


1 


1 


The  radial  location  for  zero  shear  was  used  rather  thin  the  more  frequently 
used  radius  for  maximum  velocity  because  it  is  more  sharply  defined  and  also 


applies  to  the j tangential  case  where  the  correlations  used  are 

i  -  \ 


'r0 


wi  (rxr0=  0  *  ^ 


or 


'rP 


W0  (r0  -  rT  =  05 
s  r0  u 


(36) 


The  jtwo  abscissas  are  the  fractional  distances  from  either  wall  ,:o  the  zero 

shear  location.  In  effect,  the  inner  scale  is  stretched  disproportionately  since 

r  _  .  -  r j_  is  typically  0.3  (r0  -  r^).  In  short,  the  outer  profiles  are 

tiie  more  important  since  they  represent  the  major  fraction  of  the  flow.  The 

i  j 

dashed  line  in  each  figure  represents  the  simple  mixing  length  limit, 
x  =  pd2(0u/3y)i!  where  £  =  ky  and  k  =  0.4.  Since  et/v  =  0.4  y+a 

d(c  /u  Ar)/d(y/Ar)  =  0.4.  The  constant  k  was  determined  from  fully  developed 
I  z 

pipe,  channel,  and  boundary  layer  experiments.  It  is  drawn  for  reference  purposes 
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Figure  h3 .  Dimensionless  Axial  Eddy  Diffusivity — X/D^=  9.0 
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Dimensionless  Axial  Eddy  Diffusivi'ty — X/Dh=  10.3' 
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only.,  although  it  is  quite  surprising  how  close  the  slope  of  0.4  matches  much  of 
the  data- -both  axial  and  tangential- -similar  to  the  findings  of  Liu,  et  al.  [66]. 
They  observed  that  for  a  fourfold  increase  in  surface  roughness,  the  Karman 
constant  is  unaltered,  and  the  value  e^/u^fi  is  0.09  for  the  outer  half  of  a 
turbulent  boundary  layer. 

Inspection  of  the  limiting  slopes  in  Figs.  42-46  reveal  that  for  the  axial 
diffusivities  the  Karman  constant  is  unaltered  by  swirl  for  $  =  0°  -  45°.  Values 
of  e ,/u*Ar  =0.1  apply  in  the  central  core  of  the  annulus  for  swirl  angles  of 
$  =  0“  and  10°.  Larger  axial  diffusivity  values  appear  whenever  the  local  flow 
is  more -nearly  the  free-vortex  type  (wr  =  constant).  That  is,  for  low  swirl 
rates,  0°  and  10°,  the  inlet  produced  a  nearly  solid  body  rotation  which  does 
not  alter  the  axial  diffusivities.  For  $  =  30°  and  45°,  the  early  x/D^  stations 
have  a  free  vortex  character  and  relatively  large  outer  flow  diffusivities.  The 
^  =  60°  profiles  behave  differently.  The  free  vortex  character  is  suppressed 
at  the  early  stations  due  to  strong  inlet  pressure  forces  but  is  observed  at 
the  later  stations. 

The  inner  flow  tangential  diffusivity  data  depicted  in  Figs.  47-51  chal¬ 
lenge  the  concept  of  mixing  length.  Expected  positive  slopes,  derQ/dr  later 
become  negative.  Spikes  in  the  erQ  curves  disappear  at  larger  x/D^.  For  the 
largest  swirl  rates,  the  diffusivities  approach  zero--a  result  in  agreement 
with  solid  body  flows.  There  is  no  clearly  defined  inner,  tangential  boundary 
layer.  The  velocity  in  this  region  is  already  increasing  with  a  near  forced- 
vortex  character. 

Spikes  also  occur  in  the  outer  flow  tangential  viscosity  data.  Large 
viscosities  are  observed  for  the  low  swirl  angles.  For  =  60°,  the  outer  core 
flow  viscosities  agree  remarkably  well  with  the  universal  value  of  erg/wo*ro=0.028 
derived  by  Kinney  from  Taylor's  free  vortex  experiments  with  rotating  cylinders 


Dimensionless  Tangential  Eddy  Difi’usivity- 


Dimensionless  Tangential  Eddy  Diffusivity- 
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(in  theabsenceofaxial  "carrier"  flow).  Near  the.  outer  .wall,  the  slopes 
decrease  sharply- with  increased  swirl,  indicating  that  the  Karinan  constant, 
k  =  0.4  is  sensitive  to.  swirl,  especially ^at  the  -  smaller  swirl  rates . 

It  appears  that  large  deviations  in  turbulent  diffusivities  occur  in 
regions' where  the  character  .of  the  tangential  velocity  .profile  is  changing 
from  an-  essentially  free  vortex  flow  to  an  essentially  forced  vortex  type  of 
flow.  The  observed  diffusivity  values  exhibit  unexpected  peaks  temporarily  as 
the  gradient  in  .angular  .momentum  decreases  from  0  to  negative  values.  This 
is  not  a  sharply  defined  quantity.  A  possible  classification  parameter  is 
discussed  in  Appendix  B. 

.  1 

The  present  method  of  determining  turbulent  transport  coefficients  involves 
integration  and  differentiation  of  large  amounts  of  experimental  profiles— ■ 
a  procedure  which  virtually  guarantees  a  large  amount  of  scatter.  Secondly, 
some  turbulence  intensity  terms  were  omitted  in  the  integral  equations  because 
no  equipment  was  available,  in  order  to  test  the  validity  'of  these  simplified 
equations',  a  conplete  static  pressure  field  for  all  six  of  the  Reynolds  stresses, 
and. a  complete  mean  velocity  field  for  all  three  components; are  all  needed. 

CONCLUSIONS 

On  the  basis  of  the  experimental  results,  the  following  remarks  can  be 
made  (for  the  swirl  rates  tested): 

1)  The  axial  velocity  profiles  for  $  =  0°  to  45°,  at  the  first  station 
had  the  characteristic  shape  of  a  fully  developed  profile  because  of 

the  entrance  conditions.  At  $  =  60°,  the  flow  is  approaching  a  reversed 
flow  condition  along  the  inner  wall  and  hence  the  maximum  velocity 
point  has  moved  out  closer  to  the  outer  wall. 

2)  For  $  =  10°  to  45°,  very  little  change  was  noted  in  the  axial  velocity 
profiles  with  swirl  relative  to  those  for  straight  flow. 
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The  major  characteristic  noted  in  the  decay  of  the  tangential  velocity 
for.'4>  =  10°  to  45®  was  the  change  from  a  free  vortex nature  in  the 
inlet  to  that  of  a  forced  vortex  nature  at  the  outlet.  Also,  the 
decay  did  not  occur  in  a  smooth,  gradual  fashion.  For  J  =  60°,  the 
forced  and  free  vortex  nature  of  the  profiles  did  not  change  much"--- 
with  axial  distance  and  the  decay  did  occur  in  a  smooth,  gradual 
fashion. 

The  flow  angle  for  $  =  10?  to  45°  close  to  the  inner  wall  exhibits 
a  "snake-like"  nature  as  one  moves  downstream.  This  was  not  noticed 
at  the  outer  wall  to  any  extent.  At  =  45°,  the  path  followed  by  the 
flow  close  to  the  inner  wall  would  be  similar  to  that  followed  if  one 
could  imagine  himself  walking  along  the  coils  of  a  spring  that  was 
in  various  degrees  of  tension  alternating  with  conpression  along  its 
length.  No  "snake-like"  nature  was  noted  in  the  flow  angle  with 
axial  distance  at  the  highest  swirl  rate  tested. 

Good  to  relatively  poor  agreement  was  obtained  in  the  inner  axial 
velocity  profile  between  the  data  and  the  universal  velocity  distri¬ 
bution.  The  f  =  60®  data  showed  the  poorest  agreement. 

No  regular  shift  in  the  local  axial  skin-friction  distribution  at  the 
inner  wall  is  noted  with  increasing  swirl  rate  for  <F  =  0°  to  45°.  A 
definite  decrease  in  the  axial  skin- friction  seems  to  occur  in  the 
inlet  region  for  ?  =  45°.  To  generalize,  no  significant  change  in 
axial  skin-friction  coefficient  was  found  with  swirl  rates  up  to 
$  =  30°.  The  $  =  60°  distribution  shows  a  definite  decrease  in  the 
axial  skin-friction  coefficient. 

It  appears  as  though  the  axial  pressure  gradient  will  approach  a  fully 
developed  state  before  the  swirl  conpleteiy  decays . 
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8)  Five  distinct  zones  were  observed'  in  the  tangential  velocity  profile. 

9)  Pitot  cylinder  and  hot  film  probe  velocity  profiles  are  in  agreement 
in  the  central  annular  core.  The  pitot  cylinder  indicates  high 
velocities  near  solid  boundaries. 

10)  Preston  tube  and  Clauser  chart  determinations  of  axial  and  tangential 
shear  conponents  agreed  within  IS  percent  for  the  inner  wall  shear. 

11)  Axial  and  tangential  shear  stress  profiles  were  computed  from  integral- 
differential  forms  of  the  conservation  equations  using  the  measured 
velocity  profiles  in  the  integrands.  Turbulence  intensity  terms  were 
not  measured  and  were  excluded. 

12)  Axial  and  tangential  diffusivities  were  conputed  from  the  appropriate 
shear  profiles  conputed  from  the  integral-differential  equations  and 
a  Boussinesq  hypothesis  relating  the  turbulent  Reynolds  stress  to  the 
rate  of  deformation.  In  general,  the  location  of  zero  shear  stress 
did  not  occur  where  the  rate  of  deformation  was  zero  and  their  ratio, 
the  diffiisivity,  approached  infinity-both  positively  and  negatively. 

13)  Axial  diffusivities  are  weakly  influenced  by  swirl.  This  is  probably 
due  to  the  fact  that  the  axial  velocity  profile  is  only  mildly  influ¬ 
enced  by  swirl. 

14)  For  axial,  diffusivities,  the  Karman  constant  is  unaltered  by  swirl. 

15)  The  tangential  diffusivities  behave  in  much  the  same  manner  as  those 
observed  in  curved  channels. 

16)  At  large  swirl  rates  where  the  tangential  swirl  has  a  free  vortex 
character,  the  outer  flow  diffusivities  agree  with  the  universal 
value  erQ/w*0r0  =  0*028  derived  by  Kinney  from  Taylor's  free  vortex 
experiments  with  rotating  cylinders. 


17)  Large  spikes  in  the  tangential  viscosities  occur  in  regions  where  the 
character  of  the  tangential  velocity  profile  is  changing  from  essen¬ 
tially  free  vortex  to  forced  vortex  nature. 

18)  Kinney's  conclusion  that  states  that  in  plane  curved  turbulent  flows 
which,  possess  universal  similarity,  uie  angular  velocity  is  a  trans¬ 
ferable  quantity  in  the  same  sense  as  linear  momentum  is  a  transferable 
quantity  in  rectilinear  flows ,  also  applies  for  nonsimilar  profiles 
with  an  axial  carrier  flow. 


APPENDIX  A.  Comparisons  Between  Pitot  Cylinfer  1  .-r.be  and  Hct  Film  Probe 

The  hot  film  probe  used  was  purchased  from  Thera.'  Systems,  Inc.,  St.  Paul, 
Minnesota.  A  h(  film  sensor  was  selected  instead  v".  a  hot  wire  sensor  because 
lie  film -sensor  is  more  rugged  and  less  susceptibl  i  to  fouling.  The  model  nuni>er 
of  the  probe  was  1270-10-6.  The  diameter  of  the  -  (in  sensor  was  0.001  inches  and 
the  sensing- length- was  0.020  -inches.  The  anemone  -ey  system  used  was  a  Theimo- 
Systems,  Inc.  constant  tenperature  anemometer,  'ioc'^l  1050.  The  instructions 
for  operating  the  anemometer  are  given  in  the  i  'St  ruction  manual  for  the  Model 
1050  anemometer. 

A  Beckman  voltage  to  frequency  converter ,  ^odel  651,  and  frequency  counter. 
Model  6148,  were  connected  to  the  output  of  the  anemometer  to  monitor  the  voltage 
readings.  The  instrument  was  accurate  to  within  0.3%.  This  setup  provided  an 
easy  to  road  digital  output  and  the  capability  of  averaging  the  input  signal, 
in  this  case  the  output  of  the  anemometer »  over  a  ten  second  gate  time. 

Typical  axial  and  tangential  velocity  data  gathered  with  the  two  probes 
are  presented  in  Figs.  52  and  53.  The  hot  film  flow  angles  could  not  be  deter¬ 
mined  as  accurately  with  the  hot  film  probe  as  they  could  be  with  the  pressure 
probe  (see  Figure  53) .  The  resolution  obtainable  with  the  hot  film  probe  was 
not  as  great  as  that  of  the  pressure  probe.  The  flow  angles  determined  with 
the  hot  film  probe  are  believed  to  be  within  ±1°  at  best  and  ±4°  at  worst.  Most 
of  the  flow  angles  determined  with  the  two  probes  agree  within  ±2°.  The  major 
ity  of  the  flow  angles  determined  with  the  pressure  probe  are  believed  to  be  • 
within  ±1/2°  or  better.  Very  close  to  the  inner  wall  at  high  swirl  rates, 
i.e.,  at  $  =  45°,  it  is  believed  that  the  pressure  probe  gives  false,  high 
readings,  In  equation  7,  it  is  shown  that  the  radial  pressure  gradient  is  dir¬ 
ectly  proportional  to  the  square  of  the  tangential  velocity  divided  by  the 
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radius 


w2  _  1  3p 
r  p  3r 


(37) 


As  a  -.result,  at  high  swirl  Tates  the  resulting  downwash  or  crossflow  along  the 
probe  and  the  nature  of  the  flow  around  the  probe  is  probably  quite  different  than 
that  occuring  with  weak  swirl.  Near  the  inner,  wall,  the  smaller-dimensioned 
hot  film  probe  should  not  disturb  the  flow  as  severely  as  the  pressure  probe  or 
be  affected  by  the  downwash.  Because  of  these  reasons,  it  is  felt  that  the 
pressure  probe  gives  false  readings  and  that  the  hot  film  values  are  closer  to 
the  actual  values  in  this  instance. 

The  total  velocity  measured  with  the  pressure  probe  near  the  inner  wall  at 
^  =  45°  was  also  somewhat  higher  than  that  measured  with  the  hot  film  probe. 
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APPENDIX  B.  Approximate  Locus  of  Maxinam  Turbulent  Diffusivity 

If  the  flew  can  be  written  in  the'  form  w  =  ci®,  than  r/u*  du/dr  =  m,  where 
m  =  +1  for  forced-  flew  and  m  =  -1  for  free  flow;  we  note  that  m.=~  0  (dw/dr  -  0) 
is  a  useful  boundary  between  the  two  flows.  Ihe  location  of  zero  swirl  velocity 
gradient  is  not,  however,  the  location  where  viscosity  spikes  are  observed. 

Li  classifying  this  character,  neither  the  rate  of  deformation  or  the 
vorticity  gives  a  sharp  change.  Consider  the  followirg  table 


hfixed  Flow  Classification 


Flow  Type 

Free  Vortex 

Forced  Vortex 

Velocity  w  = 

Cx/r  >  0 

C2r  >  0 

Velocity  Gradient  dw/dr  = 

-C/r2  <  0 

c2  >  0 

Deformation  Rate 

D  =  dw/dr  -  w/r  = 

-2^/r2 

0 

Vorticity 

F  =  dw/dr  +  w/r  = 

0 

2C  >  0 

2 

Product  of  D*F 

DF  =  (dw/dr  -  w/r)  (dw/dr  +  w/r) 

=  [ (dw/dr) 2  -  (w/r)2] 

<  0 

>  0 

Inspection  of  the  present  turbulent  diffusivity  data  reveals  that  the  local 
maximums  occur  close  to  the  position  where  the  product  DF  is  equal  to  zero.  The 
diffusivity  data  appear  to  be  most  accurate  outside  these  two  points  -  closer  to 


the  wallf. 
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